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ABSTRACT 


This  report  summarizes  the  work  conducted  at  the  Fhilco  Research 
Laboratories  on  Contract  AF  04(611)-9705 ,  "An  Investigation  and  Demonstration 
Of  Film  Protected-Convectively  Cooled  Nozzles".  Detail  analytical  and  experi¬ 
mental  studies  were  conducted  during  this  reporting  period  on  internal  convec¬ 
tion  and  conduction  heat  transfer,  coolant  selection,  thermochemistry,  and  the 
fluid  mechanic  aspects  of  injection  through  discrete  holes  in  a  reacting  turbu¬ 
lent  boundary  layer.  Laboratory  experiments  were  conducted  to  provide  insight, 
data,  and  direction  to  the  analysis. 
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SECTION  1 


INTRODUCTION 


The  emphasis  In  present  and  future  propellants  for  solid  rocket  motors 
has  been  toward  higher  Impulse.  This  can  be  accomplished  with  either  reducing 
the  molecular  weight  of  the  exhaust  products,  Increasing  theQflame  temperature, 
or  both.  Presently,  flame  temperatures  of  the  order  of  6500  F  to  6800°F  are 
being  considered,  and  Indeed  many  tests  have  been  conducted  with  propellants  in 
this  category.  However,  this  increase  in  temperature  has  produced  problems  in 
nozzle  design  and  integrity,  especially  for  long  duration  firings. 


The  investigations  and  results  of  contract  AF  04(611)-8387 ,  "Applied 
Research  for  Advanced  Cooled  Nozzles"  showed  that  several  thermal  protection 
concepts  are  capable  of  containing  these  high  energy  corrosive  propellants.  One 
of  these  concepts,  film  protected-convectively  cooled  nozzles,  looked  extremely 
promising,  however,  further  analytical  definition  was  required  as  well  as  experi¬ 
mental  verification. 


This  report  summarizes  the  work  conducted  during  the  first  quarter  at 
the  Philco  Research  Laboratories  on  the  "Investigation  and  Development  of  Film 
Protec ted-Convectively  Cooled  Nozzles".  This  effort  includes  an  analytical  and 
experimental  investigation  of  the  fluid  mechanic  aspects  of  injection  into  tur¬ 
bulent  boundary  layers  associated  with  reacting  nozzle  flows;  a  thermochemical 
analysis  of  the  reactions  encountered  in  the  boundary  layer  between  the  injection 
fluid  and  free  stream;  a  thermochemical  analysis  and  experimental  study  of  the 
reactions  between  the  injectant  fluid  and  the  nozzle  wall  material,  as  a  function 
of  wall  temperature.  Heat  transfer  analysis  for  both  the  internal  convection  and 
conduction  were  conducted  for  various  coolants,  heat  sink  thicknesses,  coolant 
flow  rates  and  coolant  hole  size  and  spacing.  A  very  versatile  heat  transfer 
program,  existing  at  the  Philco  Research  Laboratories,  was  used  for  the  internal 
convection  and  conduction  studies.  The  effects  of  dissociation  and  particle 
flow  with  radiation  were  also  studied  analytically  for  the  coolant  flow.  The 
flame  side  forced  convection  heat  transfer  to  the  nozzle  wall  without  film  pro¬ 
tective  was  also  investigated.  Future  investigations  will  include  the  effects 
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of  boundary  layer  injection. 

Also  during  this  quarter,  laboratory  experiments  were  conducted  in 
order  to  provide  insight  into  the  thcrmochemical  aspects  of  fluid  injection 
and  convection.  Cold  flow  experiments  for  fluid  injection  into  turbulent  bound¬ 
ary  layers  were  also  performed  in  order  to  gain  a  fundamental  insight  into  the 
phenomena  occurring,  and  Lo  provide  some  limited  data  on  flow  separation  para¬ 
meters,  mixing  region  length,  hole  size  and  spacing  variables,  as  well  as  hole 
angles  and  flow  rates. 

The  culmination  of  the  various  aspects  of  this  investigation  is  to 
provide  a  sufficient  analytical  and  experimental  background  to  allow  one  to 
reliably  design  an  optimized  film  protected-convectively  cooled  nozzle.  How¬ 
ever,  as  shown  in  this  report,  a  great  deal  of  effort  is  still  required  in 
several  of  the  above  mentioned  areas.  These  will  be  more  fully  explored  during 
the  course  of  this  program. 
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SECTION 


2 


STUDY  AND  ANALYSIS 


2.1  THERMOCHEMISTRY 


During  the  first  quarter,  a  considerable  effort  has  been  devoted  to  the 
thermochemical  evaluation  of  potential  injection  fluids.  A  preliminary  cooling 
system  design  and  weight  analysis  was  performed  to  determine  the  relative  ad* 
vantages  of  working  fluids  which,  in  their  storage  state,  would  be  gases,  liquids 
or  solids.  (This  analysis  will  be  refined  and  presented  in  a  later  report.)  It 
was  concluded  that  the  gas  storage  system  would  be  the  most  practical  for  the 
present  purposes  and  that  it  is  competitive  with  liquid  and  solid  storage  systems 
on  a  weight  basis.  Several  liquids,  particularly  NH.  and  CC1,,  apparently  have 
significant  potential  but  no  attractive  solids  were  round.  In  general,  the 
gases  and  liquids  of  interest  have  been  examined  according  to  the  following 
thermochemical  criteria: 

1.  Thermal  stability  of  gases  under  the  predicted  pressures  and 
temperatures  within  the  convection  cooling  passages 

2.  Chemical  stability  of  gases  and  their  dissociation  products 
relative  to  the  pyrolytic  .raphite  walls  of  the  convection 
cooling  passages 

3.  Chemical  stability  of  the  pure  injection  gases  relative  to  the 
pyrolytic  graphite  throat  surface  at  the  prevailing  nozzle 
pressures  and  over  the  graphite  surface  temperatures  of  Interest 

4.  Chemical  stability  of  the  injection  gases  relative  to  the 
propellant  combustion  products  as  a  function  of  pressure, 
temperature  and  mixture  ratio 
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5.  Chemical  stability  of  the  injection  gas-combustion  product  mixture 
relative  to  the  pyrolytic  graphite  throat  contour. 

In  addition,  the  following  specific  considerations  are  regarded  as 
desirable  and  have  provided  basis  for  selection  or  rejection  of  some  cooling 
fluids : 

1.  Only  endothermic  dissociation  and  no  reactions  with  the  graphite 
in  the  cooling  passages  are  allowed. 

2.  The  injection  gas  must  be  less  corrosive  than  the  propellant  com¬ 
bustion  products  at  the  graphite  design  surface  temperatures. 

3.  It  is  desirable  to  add  unstable  carbon  based  molecules  to  the 
boundary  layer. 

4.  The  injection  gas  should  have  high  potential  as  a  convection 
coolant  and  a  high  total  enthalpy  change  from  ambient  to  flame 
temperatures. 

5.  Readjustment  of  the  combustion  product  boundary  layer,  due  to 
dilution  and  consequent  chemical  reaction,  should  not  be  exother¬ 
mic  . 

Potential  Injection  fluid  classes  were  examined  and  the  following  com¬ 
pounds  were  chosen  for  detailed  study:  He,  Ar,  N.,  CO,  H-,  CH^,  BF. ,  CC1,  ,  CF^ 
and  CS-.  The  helium^argon  will  have  the  same  influence  on  nozzle  corrosion, 
with  tneir  potential  (1)  as  a  convection  coolant,  (2)  Influence  on  chemical  diffu¬ 
sion,  and  (3)  storage  weight  depend  on  their  molecular  weights.  The  nitrogen  and 
carbon  monoxide  represent  inert  constlutents  of  the  combustion  products  which 
will  behave  like  the  He  and  A  except  for  their  influence  on  the  boundary  layer 
equilibrium  composition.  Nitrogen  has  been  examined  in  detail  and  will  be  dis¬ 
cussed  separately.  The  net  thermal  effect  expected  from  any  of  the  Inerts,  when  the 
they  dilute  the  combustion  product  boundary  layer  gases,  is  slightly  endothermic 
due  to  an  increase  in  the  dissociation  of  hydrogen.  The  inerts  achieve  corrosion 
control  primarily  through  diffusion  boundary  layer  thickness  and  conductance  changes. 
The  corrosion  control  and  internal  convection  requirements  apparently  conflict  for 
the  inerts.  However,  helium  is  favored  among  these  because  of  lower  storage  weight 
and  high  impulse,  at  least  for  propellants  which  yield  the  smallest  amounts  of  CO. 
and  H20.  2 


Carbon-sulfur  compounds,  such  as  carb  n  disulfide,  CS.,  were  considered 
but  were  rejected  on  the  basis  of  the  relatively  high  decomposition  pressures  of 
carbon  monosulfide  and  gaseous  sulfur,  that  is,  2CS2 (g)"2CS(g)+S2 (g) .  At  the 
temperatures  of  interest  the  product  gaseous  sulfur  should  react  extensively  with 
the  graphite  walls  of  the  injection  capillaries  to  form  additional  gaseous  carbon 
sulfides.  The  same  conclusion  was  extended  to  cover  all  sulfur  compounds  and  this 
class  has  been  eliminated  from  further  consideration. 
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Carbon  tetraf luoride  was  considered  but  was  also  rejected  on  the  basis 
of  relatively  high  decomposition  pressures  from  the  reaction  CF^(g)»CF  (g)+F(g). 

The  monatomic  fluorine  should  react  extensively  with  the  capillary  walls  to  form 
gaseous  C^F., .  Because  the  CF^  is  the  most  stable  carbon-fluorine  compound,  the 
freons  can  also  be  eliminated.  However,  BF  because  of  its  high  thermal  stability, 
was  investigated  further.  It  was  subsequently  concluded  that  the  addition  of 
fluorine  compounds  invariably  leads  to  exothermic  reaction  in  the  combustion  pro¬ 
duct  boundary  layer.  Furthermore,  these  reactions  produce  species,  such  as  BOF, 
which  will,  in  turn,  react  with  the  graphite  walls.  11F  is  the  least  reactive  of 
the  class  and  would  be  expected  to  behave  like  an  inert  gas  of  intermediate 
molecular  weight. 

Carbon  tetrachloride  was  also  considered  since  it  decomposed  completely 
to  carbon  and  chlorine.  The  resultant  chlorine  is  relatively  unreactive  with 
the  walls  of  the  injection  capillaries.  Carbon  tetrachloride,  however,  is  a 
liquid  under  ambient  conditions;  and  its  use  involves  the  problems  of  flow 
regulation  and  boiling.  Because  of  this,  interest  in  the  use  of  CCl^  as  an 
injection  fluid  has  been  deferred  until  the  study  of  other  promising  fluids, 
which  are  and  can  be  stored  as  gases,  is  completed.  Compounds  such  as  chloro¬ 
form  (CHCl^)  are  also  of  interest.  These  will  be  investigated  further  during  the 
second  quarter  and  their  potential  advantages  will  be  pointed  out  in  the  discussion 
of  methane.  It  should  be  noted  that  the  saturated  carbon-chlorine,  carbon-hydrogen- 
chlorine  and  carbon-hydrogen  compounds  all  undergo  endothermic  dissociation,  pro¬ 
ducing  condensed  carbon.  The  saturated  compounds  tend  to  be  exothermic.  The  more 
complex  saturated  liquid  C-H-Cl  and  C-H  compounds  are  of  interest  but  their  behavior 

can  be  inferred  from  that  of  CC1,  and  CH. . 

4  4 

The  results  of  the  detailed  calculations  for  methane  and  nitrogen  are  pre¬ 
sented  in  the  subsequent  paragraphs.  The  behavior  of  hydrogen  can  be  inferred  from 
that  of  the  dissociated  methane.  Hydrogen  has  been  rejected  because  of  its  high 
reactivity  with  graphite  at  temperatures  above  about  4000  F.  On  the  other  hand, 
methane  appears  to  be  the  most  attractive  injection  gas  discovered  to  date.  This 
is  basically  due  to  the  addition  of  excess  carbon  to  the  boundary  layer.  Both  the 
combustion^product  oxidizers  and  the  hydrogen  (which  forms  C^H.  at  temperatures 
above  4000  F)  react  with  this  carbon  instead  of  the  wall  material.  Endothermic 
reaction  enhances  thermal  control  and  it  becomes  possible  to  virtually  eliminate 
contour  corrosion. 

a.  Methane 


Figure  1  shows  the  computed  equilibrium  compositlgn  of  methane  at  an 
assumed  injection  pressure  of  700  psia.  Above  1600  K  (2880  R)  the  methane  is 
essentially  all  decomposed.  Condensed  carbon,  which  is  not  Included  in  the  graph, 
is  the  major  product  of  decomposition.  As  discussed  in  Section  3.1,  Laboratory 
Experiments,  the  particle  size  of  the  condensed  carbon  has  been  found  to  be  on  the 
order  of  0.1  micron  and  should  not  significantly  lag  the  gas  flow.  To  a  large  ex¬ 
tent,  then,  the  product  hydrogen  in  the  methane  stream  remains  saturated  with  respect 
to  these  carbon  particles,  preventing  reaction  with  the  capillary  walls. 

The  point  of  injection  is  assumed  to  be  upstream  of  the  throat  at  a 
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TEMPERATURE  (DECREES  RAW  I  NT) 


FIGURE  1.  EQUILIBRIUM  COMPOSITION  OF  CHA  AT  700  PSIA  (GASEOUS  SPECIES) 
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pressure  approx imal ing  -.01)  psia  because  the  injected  gas  film  has  a  chemical 
composition  which  differs  from  that  of  the  combustion  free  stream,  It  will  act 
as  a  sink  for  chemical  species  As  the  two  streams  flow  away  from  the  point  of 
injection  they  will  tend  to  diffuse  into  each  other  because  of  the  gradients  In  the 
concentration  of  the  various  gaseous  species.  The  depth  of  penetration  of  the 
reactant  gases  into  the  injected  film  will  be  an  exponential  function  of  time, 
mass  diffusivitv,  and  any  initial  difference  in  concentration.  To  simulate  this 
variable  condition  three  ratios  of  injection  gas  to  combustion  gas  mass  flows 
have  been  considered,  1  9,  1  c,  and  91  The  combustion  product  spectrum  chosen 
is  that  from  the  6720  F  aluminized  propellant  system  simulated  in  composition  by 
the  so-called  "Hot-Solid"  composition  given  in  Table  I.  Table  I  also  lists  the 
important  products  to  be  found  in  the  free  stream  at  the  throat. 

As  evidenced  in  Figures  2  through  8,  the  equilibration  of  methane  with 
the  free  stream  sut  *tant  ia;  1  y  r»  duces  the  c  cn'f  ent  rat  ion  of  the  principal  species 
which  can  react  with  the  nozzle  wall  v.O^  •  H  O,  0,  OH,  and  condensed  A1  0  ,  with 
consequent  increases  in  the  CO  H.  H.  and  condensed  AN  species.  Figur^  9  shows 
the  disproportionate  increase  in  the  number  of  gaseous  moles  per  100  grams  of  the 
mixture  which  is  the  result  of  hydrogen  addition  and  chemical  reaction.  This 
suggests  that  there  will  be  a  substantial  amount  of  gaseous  blowing  in  the  bound¬ 
ary  layer.  In  Figure  10,  the  amount  of  wall  material  required  to  theoretically 
saturate  one  gram  of  the  gas  mixture  (injection  gas  plus  combustion  products)  is 
designated  "B".  A  negative  "B"  value  implies  that  there  is  more  than  sufficient 
carbon  from  the  decomposition  of  the  injectant  methane  gas  to  saturate  the  mixture. 
Thus,  below  3600  K  (6-»SG  R  all  methane  mixtures  should  substantially  reduce  Ihe 
extent  of  chemical  reaction  of  the  nozzle  wall  with  the  combustion  product  gases. 
The  rise  in  "B"  values  above  aoout  2500  K  (4500°R)  which  occurs  with  the  addition 
of  methane  can  be  attributed  to  the  increasing  stability  with  temperature  of 
gaseous  Until  the  "B"  values  become  positive,  at  higher  temperatures,  the 

hydrogen  resulting  from  the  CH  decomposition  will  consume  the  excess  carbon  in 
the  boundary  layer  endothermically  but  otherwise  to  no  advantage.  The  carbon- 
chlorine  injection  gases  should  lead  to  negative  "B"  values  to  significantly 
higher  graphite  temperatures,  Evidently  a  30%  CH  mixture  will  still  be  inert 
to  the  graphite  wall  at  5000°F  * 

"B"  values  have  also  oeen  computed  for  the  case  of  the  methane  within 
the  injection  capillaries,  a  d  are  shown  in  Figure  11.  Negative  "B"  values  are 
encountered  up  to  3500  K  (6270  R;  Above  this  temperature,  the  stability  of 
C2H2  gives  rise  to  positive  "E"  values. 

The  enthalpy,  in  B T / 1  b ,  of  the  free  stream  and  the  free  stream  plus 
various  weight  percentages  of  methane  have  also  been  computed  as  a  function  of 
temperature  and  are  shown  in  Figure  12  The  enthalpy  values  shown  represent  the 
summation  of  the  heat  of  formation  at  536  R  plus  the  change  in  enthalpy  above 
that  temperature  of  all  of  the  combustion  products  and  of  the  products  of  the 
reaction  of  the  combustion  products  with  methane.  These  data,  then,  Include 
the  heats  of  any  reactions  which  occur  when  methane  is  mixed  with  the  free 
stream.  Figure  12  also  gives  the  entbolpy  values  if  no  chemical  reactions  be¬ 
tween  the  methane  and  the  combustion  product  species  are  allowed  to  occur. 
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FIGURE  6.  EQUILIBRIUM  CO  CONCENTRATION  FOR  "HOT  SOLID"  -  CH4  MIXTURES 
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TABLE  I 


"HOT -SOLID"  PROPELLANT 
CALCULATED  COMPOSITION  OF  CHAMBER  AND 
EXHAUST  PRODUCTS 
(mols/100  grams  propellant) 
(PARTIAL  LIST) 


Spec ies 

Chamber 

Throat 

A1C 1 

0.0536 

0.0445 

CHO 

0.0008 

0.0006 

CO 

1.2623 

1.2632 

C°2 

0.0757 

0.0751 

Cl 

0.0568 

0.0568 

HC1 

0.1357 

0.1487 

CIO 

0.0001 

0.0000 

H 

0.2586 

0.2412 

HN 

0.0002 

0.0001 

OH 

0.0811 

0.0674 

H2 

0.4094 

0.4209 

rc 

ro 

O 

0.2090 

0.2086 

N 

0.0003 

0.0002 

NO 

0.0163 

0.0123 

N2 

0.5340 

0.5362 

0 

0.0282 

0.0215 

°2 

0.0056 

0.0042 

A1203/c 


Chamber  Pressure: 

700  psla 

Throat  Pressure: 

405  psla 

(equilibrium) 

T  chamber : 

3991°K  (7180°R) 

T  throat: 

3800°K  (6830°R) 

(equilibrium) 
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Assuming  the  enthlpy  of  mixing  is  zero,  the  difference  in  the  enthalpy  values 
at  any  given  temperature  with  and  without  reaction  represents  the  net  heat  of 
reaction  at  that  temperature.  The  reactions  are  highly  endothermic  and  are  a 
function  of  the  degree  of  mixing,  pressure  and  temperature.  As  an  example,  3 
pounds  of  methane  mixed  in  equal  proportions  with  the  combustion  products  at 
5000°F  requires  the  addition  of  about  6000  3TU  to  maintain  the  temperature.  For 
the  107«  methane  mixture,  24,000  BTU  are  absorbed.  This  effect  will  significantly 
reduce  the  nozzle  wall  heat  transfer  and  is  vastly  superior  to  the  ordinary  film 
cooling  effects  because  of  the  dependency  on  the  degree  of  mixing. 

b.  Boron  Tri fluoride 


Boron  trifluoride  appears  to  be  quite  stable  with  respect  to  dissoc¬ 
iation  at  temperatures  up  to  about  2800  K  (5040  R)  and  as  such  would  not  be  expected 
to  react  with  the  graphite  capillary  walls.  This  is  confirmed  by  the  "B"  value 
calculations  (not  shown).  As  in  the  case  of  methane,  three  ratios  of  injection 
gas  to  "Hot-Solid"  combustion  gas  mixture  ratios  have  been  considered,  1:9,  1:1, 
and  9:1.  The  results  are  not  presented  in  detail  since  fluorides  have  been  re¬ 
jected  as  injection  fluids. 

The  equilibration  of  BF^  with  the  free  stream  principally  effects  the 
H^)  and  the  aluminum  oxide  species,  both  gaseous  and  condensed.  A  reduction  in 
tne  concentration  of  these  species  is  observed  with  consequent  increase  in  con¬ 
centration  of  g  seous  HBO.  BHO  ,  B01  (large),  A10F,  A1F,  A1F. ,  A1F_,  and  monatomic 
fluorine.  1  Z  J 

The  amount  of  graphite  nozzle  wall  material  required  to  theoretically 
saturate  one  gram  of  three  stream  (injection  gas  plus  combustion  products)  is 
shown  in  Figure  13.  Interestingly,  above  3200°K  the  "B"  values  for  BF^  +  com¬ 
bustion  products  are  significantly  lower  then  the  "BH  values  for  CH^  +  combustion 
products,  indicating,  from  a  thermochemical  stantpoint,  correspondingly  lower  a- 
mounts  of  reaction  with  the  nozzle  walls.  However,  the  "B"  values  can  never  be 
negative  so  corrosion  can  only  be  reduced,  provided  that  the  mass  transfer  rates 
of  the  reactive  fluoride  species  are  not  significantly  greater  than  for  the  com¬ 
bustion  products  without  injection.  Close  examination  of  Figure  14  shows  that 
the  reactions  are  exothermic  for  t’’e  temperature  range  of  interest. 

c .  Nitrogen 


Nitrogen,  of  course,  represents  an  inert  gas  and  would  not  be  expected 
to  react  with  the  graphite  walls  of  the  infection  capillaries .  The  increasing 
stability  of  gaseous  CN  above  2600  K  (4680  R),  however,  indicates  some  minor 
reaction  will  occur  above  this  temperature 

As  in  the  case  of  methane  and  BF^ ,  three  ratios  of  nitrogen  injection 
gas  to  "Hot-Solid"  combustion  gas  mass  flows  have  been  considered,  1:9,  1:1, 

9:1.  In  general,  the  equilibration  of  nitrogen  with  the  free  stream  does  not 
have  major  effects  (other  than  dilution)  on  the  concentration  of  the  combustion 
product  species.  Exceptions  are  gaseous  N^ ,  N,  NO,  H,  and  0.  The  dilution 
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tends  to  increase  the  dissociation  of  ,  H_0,  0  ,  and  CCL,  but  the  heat  absorbed 
is  apparently  much  les*-  than  for  the  case  or  methane  addition. 

Figure  15  shows  the  amount  of  wall  material  required  to  theoretically 
saturate  one  gram  of  free  stream  (injection  gas  plus  free  stream).  The  difference 
in  "B"  values  in  going  from  0  w/o  to  90  w/o  at  any  temperature  reflects 

principally  a  dilution  effect  rather  than  changes  in  the  reactivity  of  the 
gaseous  mixture  with  the  graphite 

Figure  16  gives  the  enthalpy,  in  5TU/lb,  of  the  combustion  product  free 
stream  and  the  free  stream  plus  various  ...  t  percentages  of  •  Apparently, 
the  outer  portions  of  the  bougdary  layer  will  be  effected  by  the  endothermic 
reactions,  but  a  wall  at  5000  F  would  not  be  equally  affected  due  to  the  exoth¬ 
ermic  recombination  reactions  which  would  occur 

d .  Cone lus ions 

It  is  clear  that  the  equilibrium  thermochemical  analysis  does  not 
provide  a  completely  valid  description  of  the  boundary  layer  chemistry  in  the 
rocket  nozzle.  The  assumption  of  equilibrium  is  apparently  unavoidable  but 
there  are  at  least  two  reasons  why  it  can  be  doubted.  Thus,  it  is  difficult 
to  imagine  that  the  Al^O^  particles  could  be  in  equilibrium  in  the  injection 
boundary  layer,  either  relative  to  the  condensed  carbon  for  methane  injection 
or  to  the  combustion  product-injection  gas  mixture.  In  fact,  the  concentration 
of  Al^O  in  the  boundary  layer  may  be  significantly  different  than  in  the  main 
combustion  product  flow,  whereas  uniform  distribution  has  been  assumed.  Conse¬ 
quently,  calculations  were  performed  for  methane  injection  with  the  condensed 
Al^O^  removed.  The  results  arc-  not  completely  reduced  for  presentation  in  this 
report.  Apparently  the  removal  of  the  A^O^  ^as  little  influence  on  effective¬ 
ness  of  the  methane.  The  H^O.  CO^ ,  0.,  0,  and  OH  concentrations  are  reduced  even 

farther  and  the  zero  "B"  values  are  shifted  to  slightly  higher  temperatures  at 

comparable  mixture  ratios 

It  is  also  rather  doubtful  that  the  carbon  particles  which  result  from 

the  methane  dissociation  in  the  cooling  passages  will  achieve  equilibrium  with 

the  gas  phase.  Due  to  flow  turning,  flow  acceleration,  and  boundary  layer 

turbulence  the  carbon  particles  may  not  be  distributed  uniformly  across  the 

boundary  layer.  It  would  be  preferable  to  force  the  condensed  carbon  to  form 

after  injection  into  the  boundary  layer  to  achieve  the  smallest  particle  size. 

In  general,  this  would  require  lower  injection  temperatures  or  higher  injection 

rates,  neither  of  which  is  desirable,  Evidently,  the  efficiency  of  the  injected 

carbon  as  a  chemical  reaction  filter  can  only  be  determined  by  actual  motor  firings. 

However,  it  has  been  estimated  that  the  injection  design  flow  rate  of  one  half 

percent  of  the  subscale  combustion  product  flow  yields  carbon  equivalent  to  the 

corrosion  rate  of  about  5  mils  per  second  over  40  in  of  nozzle  surface  at 

o 

temperatures  below  about  5000  F.  Since  the  no  injection  corrosion  rate  should 
not  be  more  than  about  one  mil  per  second,  it  is  concluded  that  the  methane 
injection  provides  excess  carbon  which  will  improve  the  apparent  effectiveness 
of  the  injection. 
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FIGURE  15.  "B"  VALUES  FOR  GRAPHITE  IN  EQ 

MIXTURES  AT  400  PS  I A 


o*  n2 

107.  N2 


SOLItf'  -  N2 


bOX  N2 


06633 


DE( 

3400 


4500 


90*  N2 


DEGREES 


2400 


2800 


Another  point  whkh  should  be  considered  involves  the  assumed  compos¬ 
ition  of  the  combustion  products  which  are  displaced  by  the  injection.  Since 
solid  propellant  motors  and  the  Philco  simulator  use  ablative  chamber  liners, 
there  will  be  a  significant  amount  of  boundary  layer  mass  addition  up  stream 
of  the  nozzle  inlet  It  is  expected  that,  for  the  present  propellant  simulation, 
there  may  be  significant  .  N  ,  and  possibly  H^O  enrichment  of  the  nozzle 
boundary  layer.  There  may  also  be  some  carbon  enrichment  resulting  from  the 
injection  of  the  phenolic  pyrolysis  products  through  the  porous  chamber  liner. 

It  is  interesting  to  spe*.ulait  that  the  corrosion  rates  with  no  injection  in 
solid  propellant  motors  may  be  significantly  effected  by  the  choice  of  ablative 
liner  material,  engine  size,  the  time  dependence  of  the  pyrolysis  rate  and  the 
use  of  porous  graphite  liners.  These  factors  should  all  be  considered.  It  is 
currently  felt  that  the  pyrolysic  mass  addition  generally  increases  the  require¬ 
ments  for  nozzle  surface  temperature  control  to  minimize  the  influence  of  the 
hydrogen  enrichment.  At  the  same  time,  it  is  expected  that  temperature  control 
will  be  potentially  more  effective  in  corrosion  control. 

In  spite  of  these  deviations  from  the  ideal  conditions  assumed  for 
analysis,  it  is  felt  that  the  control  of  corrosion  by  boundary  layer  injection 
will  be  extremely  effective  The  major  theoretical  advantages  of  carbon  inject¬ 
ion  provide  the  basis  for  preferentially  choosing  injection  fluids  from  the 
C-H,  C-Cl  or  C-H-Cl  compounds  A  number  of  inert  gases  can  provide  significant 
corrosion  reduction  through  the  chemical  sink  and  dilution  mechanisms. 

During  the  second  quarter,  work  will  proceed  in  evaluating  (1)  the 
methane  -  no  Al^O^  (2)  the  CC1,  and  (3)  the  C-H-Cl  injection  systems.  It  is 
anticipated  that  some  additionaf  calculations  will  be  required  to  provide  the 
temperature -pressure -enthalpy-compos  it  ion  relations  for  the  injection  fluid  or 
fluids  chosen  for  the  subscale  tests.  This  information  will  eventually  be 
required  for  use  with  the  boundary  layer  film  analysis. 


) 
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2.2  FLUID  MECHANICS 


A  primary  advantage  of  the  nozzle  concept  under  study  is  that  fluid 
injection  over  the  throat  contour  provides  the-  designer  with  the  opportunity  to 
moderate  the  interaction  of  the  combustion  produc t s wi th  the  nozzle  wall.  It  is 
clear  that  materials,  chemical,  thermal  and  structural  responses  will  be  effected 
and  can  be  controlled  to  the  same  degree  that  the  injection  fluid  behavior  can  be 
predicted  and  controlled.  The  boundary  layer  fluid  mechanics  studies  have  es¬ 
sentially  been  split  into  two  areas,  namely,  (1)  fluid  injection  and  (2)  injection 
film  analysis.  These  areas  will  be  discussed  briefly  before  proceeding  with  the 
presentation  of  the  results  of  the  work  of  the  first  quarter. 

A  great  deal  of  study  and  analysis  of  slot  and  hole  injection  for  *  il  i 
cooling  has  preceeded  this  investigation.  Unfortunately,  the  results  have  lacked 
sophistication  and/or  feature  mass  injection  rates  which  are  greater  than  either 
those  considered  here  or  the  displaced  boundary  layer  mass  flow.  Because  the 
present  injection  is  essentially  a  boundary  layer  problem,  very  little  of  the 
previous  work  in  film  cooling  is  thought  to  apply.  It  is  clear,  however,  that  the 
most  efficient  film  injection  can  be  achieved  by  tangential  slot  injection.  For 
the  present  design  concept,  very  low  mass  injection  rates  lead  to  very  narrow 
slots,  the  dimensional  stability  of  which  could  not  be  maintained  during  a  nozzle 
test.  A  transpiration  injection,  such  as  it  described  in  Reference  1  ,  could 
also  efficiently  generate  a  uniform  film.  This  approach  would  be  unsatisfactory 
for  coolants  which  decompose  to  carbon.  Other  disadvantages  of  this  injection 
method  were  considered  during  contract  AF  04(611) -8387  and,  consequently,  it  is 
not  actively  being  considered.  The  last  choice  is  to  inject  the  fluid  through  a 
large  number  of  discrete  holes  in  the  wall.  The  advantages  of  this  type  of  in¬ 
jection  are  currently  thought  to  outweigh  the  disadvantages.  However,  the  three- 
dimensional  wall  jet,  with  heat  transfer,  diffusion  and  chemical  reaction,  cannot 
be  reasonably  analyzed.  Therefore,  the  study  of  hole  type  injection  is  primarily 
being  approached  experimentally.  Areas  of  interest  which  have  been  or  will  be 
investigated  include: 

(1)  the  influence  of  the  various  jet  and  free  stream  flow  parameters 
on  the  jet  attachment  to  or  detachment  from  the  wall, 

(2)  the  momentum  mixing  which  occurs  along  the  jet-free  stream  boundary, 

(3)  the  influence  of  injection  parameters  on  Jet  spreading  and 
coalescence , 

(4)  the  influence  of  the  initial  jet  and  free  stream  boundary  layers 
on  the  injection  interaction, 

(5)  the  Influence  of  jet  density  and  nozzle  acceleration  on  the  In¬ 
jection  and  spreading. 

From  the  results  obtained  to  date,  it  is  apparent  that  the  injection 
jets  will  attach  to  the  wall  for  the  anticipated  nozzle  design  and  flow  conditions 
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(including  angular  injection).  It  is  also  evident  that  the  diffusive  spreading 
of  the  injection  fluid  is  slow  so  that  jet  coalescence  in  short  distances  can  be 
achieved  only  through  extensive  turbulent  mixing  (not  desired)  or  by  increasing 
the  work  done  on  the  jet  by  the  free  stream.  In  the  latter  case,  the  introduction 
of  a  small  circumferential  swirl  component,  to  the  injection  velocity  vector,  may 
spread  the  jet  without  introducing  excessive  mixing  or  turbulence.  The  recognition 
of  the  overall  complexity  of  the  injection  interaction  problem  has  delayed  (and 
possibly  precludes)  the  formulation  of  definitive  statements  concerning  the  choice 
of  parameters  such  as  hole  angles,  hole  spacing .  hole  size,  axial  position  for 
injection,  fluid  density,  etc.  However,  it  has  been  tentatively  concluded  that; 
injection  in  the  flow  acceleration  region  is  not  desirable;  minimum  practical  hole 
spacings  are  desired;  injection  angles  near  15°  are  satisfactory,  the  fluid  in¬ 
jection  velocities  should  approach  the  local  free  stream  velocity;  and  minimum 
hole  size,  approaching  the  displaced  combustion  product  boundary  layer  thickness, 
is  desirable.  It  should  also  be  noted  that,  the  generation  of  a  continuous  film 
is  not  an  absolute  requirement  of  the  injection.  Of  course,  a  discontinuous  film 
will  not  be  as  efficiert,  thermallv  or  chemically,  as  a  continuous  one.  Also, 
additional  turbulence  generated  by  the  swirl  type  injection  could  significantly 
reduce  the  protection  length  of  a  continuous  film  generated  in  that  manner. 

It  has  been  necessary  to  assume  that  the  injection  could  produce  a 
continuous  film  for  the  purposes  of  analyzing  the  subsequent  behavior  of  the 
injection  fluid  in  the  nozzle.  Immediately  after  injection,  the  two  streams 
will  tend  to  exchange  mass,  momentum  and  energy  at  rates  which  should  decrease 
exponentially  with  time.  Any  turbulence,  in  addition  to  that  characteristic  of 
the  two  streams  at  injection,  generated  due  to  the  momentum  mixing  at  the  free 
stream-jet  interface  will  tend  to  accelerate  the  mass  and  energy  mixing.  That  is, 
injection  induced  turbulent  convective  mixing  will  be  more  important  than  the  mass 
diffusion  and  thermal  conduction.  Apparently,  it  is  desirable  that  the  net  transfer 
cf  momentum  between  the  streams  be  close  to  zero.  The  extent  of  the  initial 
injection  mixing,  and  consequently  the  subsequent  thermal  and  chemical  behavior 
of  the  modified  film,  cannot  be  accurately  predicted  because  of  poor  understanding 
of  the  fundamentals  of  turbulence  Following  the  initial  mixing,  the  momentum  and 
energy  discontinuities  in  the  bourdary  layer  will  tend  to  disappear  and  the  modified 
boundary  layer  will  continue  to  develop  in  a  relatively  normal  manner.  In  general; 
the  new  energy  thicknesses  will  be  larger;  the  momentum  thicknesses  are  likely  to 
be  larger,  but  could  be  smaller  for  high  injection  mass  flux;  the  two  fluids  will 
continue  to  mix  chemically  under  the  influence  of  convective  and  diffusional 
mass  transport;  and  the  detailed  film  boundary  layer  development  will  be  dependent 
on  the  chemical  reactions  which  occur  during  mixing.  Under  reasonable  injection 
conditions , the  composition  of  the  fluid  immediately  adjacent  to  the  wall  will  not 
change  (from  the  injection  fluid  composition)  for  a  distance  downstream  which  may 
be  as  great  as  ten  hole  diameters  Subsequently,  dilution  by  the  most  mobile  of 
the  free  stream  species  will  proceed  For  the  methane  injection,  the  dilution 
by  the  reactive  species  (HO  CO,  0  ,  etc  )  will  be  retarded  or  eliminated  by 
chemical  reactions  in  the  chemical  ml  sing  region  of  the  boundary  layer.  The 
compositional  readjustments  ate  espected  to  persist  throughout  the  nozzle, 
partially  due  to  the  very  short  residence  times.  The  heat  transfer  in  the 
nossle  throat  region  will  be  reduced  by  several  mechanisms,  namely: 
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(1)  the  increases  in  boundary  layer  energy  and  momentum  thicknesses 
due  to  the  Initial  energy  and  momentum  decrements  of  the  injected 
fluid, 

(2)  similar  increases  in  boundary  layer  parameters  and  heat  blocking 
due  to  endothermic  chemical  reactions, 

(3)  repression  of  free  stream  recombination  reactions  near  the  wall, 

(9)  alteration  of  the  boundary  layer  heat  transfer  properties  due  to 
composition  changes. 

It  should  be  noted  that  the  similarity  between  the  heat  and  mass  transfer 
phenomenon  suggests  that  corrosion  rates  will  also  be  reduced  in  the  same  manner. 

It  may  also  be  observed  that  the  increases  in  boundary  layer  thicknesses  will 
induce  proportional  thrust  losses  unless  the  injection  fluid  has  higher  Impulse 
that  the  combustion  products  it  replaces.  This  would  be  a  strong  argument  for  the 
use  of  methane  (which  decomposes  to  hydrogen)  or  helium  except  that  the  boundary 
layer  losses  are  small  to  begin  with  and  the  preservation  of  the  throat  contour 
precludes  other  thrust  losses. 

Investigation  of  the  fluid  mechanic  properties  of  a.  three-dimensional 
wall  Jet  is  currently  in  progress.  This  study  is  divided  into  three  avenues  of 
approach,  all  directed  toward  the  achievement  of  an  efficient  film  injection 
design.  The  avenues  consist  of  (1)  cold  flow  and  chemical  reaction  experiments, 

(2)  mathematical  solution  of  the  pertinent  boundary  layer  equations,  and  (3)  a 
literature  survey.  Since  the  cold  flow  and  chemical  reaction  experiments  are 
discussed  in  Sections  3.2  and  3.1,  respectively,  this  section  discusses  the 
latter  two  methods  of  approach. 

During  the  first  quarter  of  the  contract  period  there  have  been  two 
mathematical  treatments  devised  to  calculate  fluid  mechanic  and  heat  transfer 
phenomena  in  the  nozzle  with  film  injection.  Each  of  these  is  regarded  as 
impractical  for  application  to  the  present  problem.  The  first  assumes  uniform 
slot  Injection  parallel  to  the  free  stream.  The  boundary  layer,  downstream  from 
injection  was  divided  into  two  regions.  In  the  region  near  the  wall  (laminar 
layer)  a  constant  shearing  stress  was  assumed  (linear  velocity  distribution)  and 
the  momentum  equation  was  solved  by  the  Integral  method.  This  layer  is  assumed 
to  be  defined  by,  and  contain,  the  injection  fluid.  In  the  second  region  (turbulent 
layer)  a  1/7  power  velocity  distribution  was  assumed  and  the  appropriate  momentum 
equation  solved.  To  combine  the  two  momentum  equations,  the  shearing  stresses 
are  equated  at  the  Interface  of  the  two  regions.  The  height  of  the  laminar  layer 
is  found  from  a  mass  balance  on  the  injected  fluid  and,  by  similarity,  the  heat 
transfer  can  be  found.  The  assumptions  of  a  linear  velocity  and  lnjectant  mass 
balance  in  the  laminar  layer  are  invalid  when  the  injected  velocity  or  the  height 
of  the  injection  slot  is  large  relative  to  the  velocity  and  thickness  of  the 
laminar  sublayer  existing  in  the  nozzle. 

A  second  analysis  was  devised  for  predicting  the  temperature  and  velocity 
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distribution  that  exists  in  the  mixing  and  developed  regions  of  f'.ie  injection 
film.  The  analysis  was  not  actually  carried  out  because  of  a  requlrem.  .t  for 
the  development  of  a  very  complicated  computer  program  for  performing  the  calculi 
tions.  A  brief  description  of  the  analysis  follows. 


Assuming  a  two-dimensional  wall  jet  and  tangential  Injection  with 
respect  to  the  surface,  the  convent  ial  mass,  momentum  and  energy  equations  for 
a  compressible,  turbulent  boundary  layer  may  be  written.  Assuming  that  the 
Lewis  and  turbulent  Prandlt  numbers  are  unity  (diffusion  and  energy  equation 
are  independent  and  eddy  viscosity  equals  eddy  conductivity),  the  boundary  layer 
is  separated  into  two  regions  and  the  eddy  viscosities  of  Dii  srler  are  used 
(Ref.  2  )  . 


Region  1  (t.aar  wall) 


Region  2 


where 
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Introducing  a  stream  function  which  satisfies 
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the  continuity  equation 
■  -  ^V)  one  may  solve 


the  momentum  and  energy  equation  by  a  numerical  method  if  the  initial  velocity 
and  temperature  profiles  are  given  at  the  point  of  injection.  It  may  be  noted 
that  no  limitation  is  imposed  on  the  velocity  distribution  in  the  jet  mixing 
region;  although  a  singularity  point  may  arise  at  the  point  of  injection,  i.e., 
where  the  free  stream  transfers  from  the  wall  to  a  free  shear  region.  This 
problem  has  been  solved  by  Denison  and  Baum  (Ref.  3)  by  use  of  a  transformation 
which  eliminates  this  singularity  point. 


At  the  present  time  a  third  mathematical  approach  to  the  injection 
film  analysis  is  in  progress.  This  method  will  eliminate  the  poor  assumptions 
of  the  first  method  and  the  excessive  cost  of  the  second. 


A  literature  searci.  ;s  in  progress  and  certain  conclusions  have  been 
obtained  concerning  various  film  parameters.  Some  of  these  conclusions  appear 
to  require  experimental  verification.  A  few  of  the  important  parameters  are 
discussed  below. 


Seban  and  Back  (Ref.  4  )  have  found  that  the  effect  of  an  accelerating 
flow  on  the  thermal  protection  aspects  of  film  injection  depends  largely  on  the 
bcation  of  large  acceleration  with  respect  to  the  initial  momentum  mixing  region. 
These  investigators  have  shown,  experimentally,  that  a  remarkable  decrease  in 
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thermal  protection  occurs  if  large  acceleration  exists  in  the  mixing  region, 
whereas  a  negligible  decrease  results  if  acceleration  is  outside  of  the  mixing 
region.  Although  the  experiments  were  performed  for  cold  flow  and  with  a  two- 
dimensional  wall  jet,  these  results  will  be  applied  to  a  test  nozzle,  in  which 
the  mixing  region  will  exist  upstream  of  an  inlet  area  ratio  of  1.3  (starting 
point  of  large  acceleration). 

The  dependence  of  the  initial  free  stream  velocity  boundary  layer  thickness 
on  heat  transfer  has  been  reported  in  Reference  5  .  It  was  experimentally  shown 
that  relatively  small  changes  in  heat  transfer  resulted  from  large  changes  in 
upstream  boundary  layer  thickness.  This  conclusion  was  applicable  to  a  two- 
dimensional  jet  and  may  not  be  applied  to  a  three-dimensional  jet.  That  is, 
because  of  the  difference  in  the  hydrodynamic  flow  field,  a  three-dimensional  wall 
jet  has  a  relatively  large  contact  area  or  region  with  the  free  stream  (in 
comparison  with  a  two-dimensional  jet)  and,  thus,  the  shearing  stress  at  the  jet- 
free  stream  interface  may  significantly  alter  the  flow  field.  Therefore,  it  is 
concluded  that  the  study  of  the  variation  in  velocity  boundaty  layer  thickness 
should  be  included  in  the  cold  flow  tests. 

The  machining  properties  of  the  nozzle  material  (pyrolytic  graphite) 
impose  certain  limitations  upon  selection  of  the  injection  hole  geometry.  These 
limitations  result  from  the  fact  that  it  is  almost  impossible  to  accurately  drill 
a  hole  less  than  0.040  inch  in  diameter  at  angles  greater  than  45°.  The  angle 
limitation  implies  that  the  injection  angle  (with  respect  to  free  stream)  may  be 
at  maximum  15°  for  a  30°  nozzle  inlet  angle. 

The  effect  of  injection  angle  on  heat  transfer  has  been  reported  in 
References  6  and  7  .  It  was  found  that  injection  angle  significantly  influences 
downstream  heat  transfer.  This  may  be  due  to  the  rather  strong  turbulent  transport 
near  the  injection  point  that  is  produced  by  angled  injection.  Apparently,  for 
maximum  thermal  protection, a  minimum  injection  angle  is  required.  Currently 
the  effect  of  injection  angle  on  wall  coverage  is  being  studied  in  the  cold  flow 
experiments  with  no  definite  conclusions  resulting  as  of  yet. 

From  the  machining  limitations,  a  hole  size  of  0.040  inch  was  selected 
for  the  nozzle.  The  selection  of  a  minimum  hole  size  is  indicated  from  the 
fact  that  for  minimum  turbulent  mixing  between  jet  and  free  stream,  it  is  desired 
that  the  wall  jet  remain  near  the  nozzle  surface  and  expose  a  minimum  amount  of  its 
surface  area  to  the  free  stream.  That  is,  minimum  turbulent  mixing  results  when 
the  jet,  essentially,  replaces  a  portion  of  the  free  stream  boundary  layer.  If 
the  jet  diameter  is  sufficiently  greater  than  the  boundary  layer  height,  the  intense 
turbulence  of  the  free  stream  will  rapdily  initiate  mixing  and  a  large  amount  of 
"protective  fluid"  will  be  lost.  This  suggests  that  for  ideal  injection  a  two- 
dimensional  slot  with  a  height  equal  to  approximately  the  laminar  sublayer  of  the 
free  stream  is  desired.  Therefore,  in  order  to  approach  slot  injection  and  reduce 
turbulent  mixing,  a  minimum  hole  size  and  spacing  are  desired.  Choice  of  hole 
spacing  is  dependent  on  jet  spreading  and  is  discussed  in  Section  3.2. 


-31- 


For  a  given  total  injectant  mass  flow  rate  and  hole  spacing,  the  sel-  tion 
of  a  minimum  injection  diameter  produces  a  maximum  (  £  U)inj  /  (  9  U)^  .  It  has 
been  found  (Ref.  5  )  that  for  a  two-dimensional  jet  the  amount  of  thermal  protection 
increases  with  (  9  U)inj  /  (  9  an<*  reaches  a  maximum  at  a  ratio  of  0.90.  Ratios 

larger  than  1.1  result  in  a  decrease  of  thermal  protection.  Therefore,  assuming 
that  this  applies  to  a  three-dimensional  jet,  it  is  desired  to  obtain  an  injection 

ratio  near  unity.  Additional  influences  of  the  U  ratio  on  jet  fluid  mechanics 
are  discussed  in  Section  3.2.  It  may  be  of  interest  to  note  that  the  9^  ratio 
for  a  single  hole  in  the  actual  nozzle  will  not  exceed  0.80  f  methane  injection. 

It  was  origionally  suggested  that  a  study  of  injecting  the  coolant  through 
a  passage  which  consisted  of  a  rece~  .  in  the  wall  would  be  of  interest.  This  design 
has  been  rejected  because  of  the  fluid  dynamic  properties  of  a  step.  For  low 
(  9  U>£nj  /  (  9  U)cd  the  fl°w  approaches  simple  step  flow.  That  is,  down  stream 
of  the  step  the  heat  transfer  increases  considerably  due  to  the  reattachment  of 
the  free  stream  to  ihe  wall.  From  Reference  3  it  has  been  shown  that  this  phenomenon 
exists  for  (  9  U)inj  /  (  9  as  ^ar8e  as  0-20  under  cold  flow  conditions. 
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2.3 


HEAT  TRANSFER 


a .  Thermal  Analyzer  Program 

Because  of  Che  three-dimensional  geometry  and  the  effects  of  the  cooling 
gas,  analytical  solutions  are  not  adequate  for  predicting  the  temperature  distri¬ 
butions  within  the  nozzle  heat  sink.  Calculations  of  temperature  distributions 
were  therefore  performed  with  the  aid  of  a  digital  computer  program  call-!  the 
Thermal  Analyzer. 

The  equations  solved  by  the  Thermal  Analyzer  can  be  obtained  by  differenc 
ing  the  general  partial  differential  equation  for  conduction  in  solids  (For  Fourier 
equation).  Conceptually,  however,  it  is  simpler  to  derive  the  equations  by  per¬ 
forming  an  energy  balance  on  a  small  section  of  the  heat  sink  material. 

The  first  step  is  to  imagine  that  "he  heat  sink  is  divided  into  a  number 
of  more  or  less  rectangular  sections.  These  sections  are  commonly  called  nodes 
since,  in  the  electric  circuit  analogous  to  a  thermal  conduction  problem,  there 
is  a  node  corresponding  to  each  section.  To  obtain  maximum  accuracy  for  a  given 
computing  time,  relatively  small  nodes  should  be  used  in  regions  where  high  tempera 
ture  gradients  are  expected.  Each  node  is  identified  by  a  number. 

The  next  step  is  to  perform  an  energy  balance  on  each  node.  The  heat 
capacity  of  a  node  times  the  node  temperature  rise  in  a  given  time  interval  must 
equal  the  heat  conducted  into  the  node  from  all  adjoining  nodes  during  tne  in¬ 
terval,  i.e. 
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number  m)  at  the  start  of  time  interval  t,°F 


R  *  the  thermal  resistance  between  nodes  m  and  n, 
m,n 

o_ 

Fsec 

BTU 

£t  ■  time  interval,  sec 


-33- 
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This  equation  can  be  rearranged  to  give 


T  1 
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-  T 
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Tm  -  Tn 

J  R 

(2) 

n 

m  m,n 

the  quantity 

m 

Tm  -  Tn 
R 

is  assumed  to  be  constant  over 

the  time 

m,n 

interval.  The  error  introduced  by  this  assumption,  called  truncation  error,  may 
be  reduced  by  using  a  small  At. 

Equation  (1)  actually  represents  a  set  of  equations,  one  for  each  node. 
The  equations  are  applied  step  by  step  starting  with  known  initial  temperatures. 
The  temperatures  computed  for  the  end  of  one  time  step  form  the  starting  tempera¬ 
tures  for  the  next  step.  Thus  complete  temperature  histories  for  each  node  are 
obtained. 


It  can  be  shown  that  the  numerical  solution  described  above  will  be 
unstable  if  the  time  increment,  £  t,  is  too  large.  The  condition  which  must  be 
met  to  ir  :  urr'  stability  is 

SeC 

n  m  R 

m,n 


for  each  node.  The  Thermal  Analyzer  Program  automatically  computes  the  "RC 

time”,  a  t  ,  for  each  node  and  sets  A  t  to  1/4  of  the  minimum  value.  This  not 
RC 

only  insures  numerical  stability,  but  tends  to  keep  the  truncation  error  small. 

The  nodal  division  used  for  one  nozzle  configuration  is  shown  in  Figure 
17.  This  configuration  will  be  used  to  illustrate  how  the  Therm*'  Analyzer  Pro¬ 
gram  is  being  applied.  Note  that,  because  of  symmetry,  the  profile  surfaces  of 
the  wedge  section  shown  can  be  assumed  insulated. 


Equation  (2)  written  for  node  4  of  Figure  17  would  be 
_  1 


T  T  T  -  T 
AJL  -  .4  +  _ _4 

4  C.  R.  ,  R-  . 

4  4  ,hg  3,4 
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(the  subscript  hg  represents  the  ho. 

The  heat  capacity  of  node  4  is  given  by 

C4  - 


-34- 


I 

u> 

Ln 

i 


I 


SCALE  (INCHES) 


R06665 


FIGURE  17.  TYPICAL  NODE  BREAKDOWN  -  SUB-SCALE  NOZZLE 


where  V,  is  the  volume  of  node  4,  (>  is  the  pyrolytic-graphite  density,  and 
c  is  the  PG  specific  heat.  Since  specific  heat  is  a  function  of  temperature, 
it  is  evaluated  at  each  time  step  as  a  function  of  T^. 

The  thermal  resistance  between  node  4  and  the  hot  gas  is  taken 


1 


4,hg 


hhgA4,hg 


where  h,  is  the  local  hot  gas  heat  transfer  coefficient  at  the  center  of  node 

4  and  A.® .  is  the  area  of  the  face  of  node  4  exposed  to  the  hot  gas. 

4,hg 

The  thermal  resistance  between  nodes  4  and  10  was  assumed  to  be 

„  .  (AA o 


4,10 


A.  i 
4,10  a 


where  is  the  distance  between  the  centers  of  nodes  4  and  10,  is 

the  area4if  the  common  interface  between  nodes  4  and  10,  and  k  is  the’PG  a- 
direction  conductivity.  Since  k  is  a  function  of  temperature,  it  was  evalua¬ 
ted  at  each  time  step  as  a  function  of  the  average  temperature  between  nodes  4 
and  10. 


The  thermal  resistance  between  nodes  3  and  4  was  computed  from 


f  3.4 
A3,4k 


c 


where  k  is  the  PG  c-direction  thermal  conductivity  evaluated  at  the  average 
temperature  between  nodes  3  and  4. 

The  formulas  given  for  C. ,  R,  ,  ,  R.  .n,  and  R~  .  illustrate  the  methods 

°  4  4,hg  4,10  3,4 

used  to  compute  all  "circuit  constants"  except  those  required  in  the  cooling  gas 
calculations. 

Because  of  the  flow  of  the  cooling  gas,  special  equations  are  required 
for  the  cooling  gas  temperature  evaluation.  (While  the  hot  gas  also  flows,  no 
special  equations  are  required  for  the  hot  gas  temperature  evaluation  since  the 
recovery  temperature  can  be  assumed  constant  with  negligible  error). 

The  cooling  gas  temperature  is  described  by  the  following  differential 

equation. 
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where  A  i9  the  local  cross-sectional  area  of  the  cooling  passage,  P  g  is  the  local 
cooling  gas  density,  c  is  the  cooling  gas  specific  heat  (function  of  temperature), 
w  is  the  flow  rate,  h  S  is  the  local  heat  transfer  coefficient,  P  is  the  local 
cooling  passage  perimefer,  Tw  is  the  local  PG  wall  temperature,  t  is  time,  and  x 
is  a  coordinate  measured  along  the  cooling  passage  axis.  An  order  of  magnitude 
analysis  shows  that,  for  the  conditions  being  investigated,  the  first  term  on 
the  left  hand  side  of  this  equation  may  be  neglected.  The  re:r... '  ..lag  equation 
may  be  differenced  to  give* 
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where 

T  is  the  temperature  of  cooling  gas  fluid  node  n,  °F 
n 

T^  +  ^  is  the  temperature  of  the  next  fluid  node  downstream  of  n,  F 
is  the  distance  between  nodes  m  and  n  (along  axis  of  cooling 
passage) 

is  the  temperature  of  a  PG  node  with  a  center  half-way  between 
nodes  n  and  n  +  1 


In  this  equation,  c  is  evaluated  at  the  average  temperature  of  nodes  n  and  n  +  1. 

g 

With  the  set  of  equations  represented  by  (4)  the  temperatures  of  all 

cooling  gas  nodes  can  be  computed  at  the  start  of  each  time  step.  The  known 
entrance  temperature  is  used  to  compute  the  first  downstream  temperature  which 
in  turn  is  used  to  compute  the  next  downstream  temperature,  etc.  In  cases  where 
the  gas  specific  heat  is  a  function  of  temperature,  it  is  necessary  to  evaluate 

c  with  the  temperature  obtained  from  the  previous  time  step  in  order  to  keep 

g 

(4)  an  explicit  equation  for  T  ,  , 
r  n  +  1 . 


*  The  cooling  passage  is  assumed  to  be  divided  into  "fluid  nodes".  In  Figure 
there  are  fluid  node  centers  coincident  with  the  center  of  node  12;  at  the  centers 
of  the  interfaces  between  nodes  7,  8,  9,  10,  11  and  12;  and  at  the  downstream  end 
of  node 
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b. 


Conduction  -  Internal  Convection 


(1)  Introduction 


This  area  of  the  analysis  is  concerned  with  the  interrelation  of 
a  pyrolytic  graphite  heat  sink  with  heat  internally  extracted  by  convection  to  a 
gaseous  coolant.  The  purpose  of  the  convection  cooling  is  to  lower  the  asymptotic 
value  of  the  wall  temperature  so  that  corrosion  will  be  decreased.  This  analysis 
will  ultimately  be  combined  with  the  injected  film  analysis  to  obtain  an  overall 
optimization.  Presently,  however,  the  study  involves  only  the  conduction  -  internal 
convection  with  given  hot  i»as  side  convection  coeilicients  to  indicate  the  heat 
sink  temperature  response  with  variations  in  heat  sink  size,  coolant  flow  rate, 
material  and  coolant  thermal  properties,  and  scaling. 

Initially,  the  problem  was  attacked  empirically  utilizing  data  obtained 
in  Contract  AF  04(611) -8387 .  Then  the  thermal  analyzer  computer  program,  des¬ 
cribed  in  Section  2. 3n ,was  applied  for  detailed  solutions.  The  purpose  of  the 
empirical  analysis  was  to  estimate  heat  absorption  required  by  the  coolant  as  a 
function  of  wall  temperature  and  duration.  This  analysis  aided  in  selecting 
inputs  for  the  computer  program. 

Since  methane  appears  to  be  the  best  coolant  for  this  application  based 

on  chemical  and  heat  transfer  studies  (See  Section  2.1  )»  determination  of  its  heat 

transfer  coefficient  in  a  tube  with  chemical  reaction,  two-phase  flow  and  rad¬ 
iation  effects  is  in  progress.  Preliminary  results  of  this  analysis  indicate  a 
significant  improvement  in  the  heat  transfer  coefficient  over  that  predicted 
using  the  properties  of  undissociated  methane.  The  computer  calculations  pres¬ 
ented  in  this  report  are  based  on  the  latter,  hence  they  are  quite  conservative. 


(2)  Empirical  Analysis  of  Cooled  Heat  Sink 


A  preliminary  approach  to  study  the  percentages  of  the  total  heat 
transferred  to  the  nozzle  which  are  absorbed  by  the  coolant  and  the  solid  heat 
sink  was  accomplished  using  the  computer  results  and  test  results  of  the  first 
demonstration  nozzle  tested  under  Contract  AF  04(611) -8387 .  These  results  are 
reported  in  Ref.  8  and  are  reproduced  here  in  Figure  18.  The  nozzle  pro¬ 

duced  5000  lb.  thrust  with  a  2.5  inch  diameter  throat  and  had  a  hgat  sink  out¬ 
side  diameter  of  7.0  inches.  The  combustion  temperature  wjsq6450  F  and  the 
convection  heat  transfer  coefficient  was  0.0053  BTU/sec  in  F.  These  vglues 
are  somewhat  lower^than  those  being  considered  in  this  contract  (T  *6720  F  and 
h  -.006  BTU/sec  in2  °F);  however,  they  are  representative  and  are  Sufficient 
f8r  tlie  purposes  of  this  preliminary  analysis. 


By  plotting  log  T  versus  log  r  with  time,  t,  as  a  parameter  where 
T  in  temperature,  °P,  and  r  is  radius,  inches;  an  empirical  equation  can  be 
written  of  the  form 


T 


a(t)  r 


b(t) 


where  a  linear  log-log  plot  is  obtained.  The  actual  plot  was  very  nearly  linear 
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FIGURE  18.  PYROLYTIC  GRAPHITE  HEAT  SINK  TEMPERATURE  HISTORY 
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hence  this  method  could  be  utilized.  The  equation  obtained  was 


T 


3550  t 


0.11 


r 

r 


12. At 


■0.65 


where  the  r  is  the  inside  radius  of  the  nozzle.  The  numerator  of  the  right 
hand  side  i§,  of  course,  the  wall  temperature.  By  holding  the  wall  temperature 
at  a  desired  value,  the  heat  sink  temperature  distribution  as  a  function  of 
time  can  be  obtained.  This  method  gives  a  result  approximately  the  same  as 
can  be  obtained  from  the  graphical  Schmidt  method.  This  method  assumes  the  sur¬ 
face  temperature  control  is  obtained  by  extracting  the  surplus  heat  at  r*r  . 

Figure  19  shows  the  result  obtained  by  allowing  the  wall  temperature  to  ?ncrease 
to  4800  F  and  then  held  constant,  as  compared  with  the  normal  heat  sink  response. 

The  heat  stored  in  the  heat  sink  at  any  time  can  easily  be  computed 
using  appropriate  temperature  profile  and  the  total  heat  versus  temperature 
curve  for  pyrolytic  graphite  shown  in  Figure  20  (Ref.  9  ).  Also  knowing 

the  wall  temperature  history,  the  total  heat  transferred  to  the  wall  up  to  any 
given  time  can  be  calculated.  The  results  ot  this  calculation  for  the  nozzle 
described  at  the  beginning  of  the  section  are  shown  in  Figure  21  .  The  heat 
which  must  be  removed  to  maintain  a  wall  temperature  of  4800  F  is  the  difference 
between  the  top  curve  (heat  stored  in  heat  sink  with  temperature  control).  The 
curve  showing  heat  transferred  to  and  absorbed  by  an  uncooled  heat  sink  is 
included  for  comparison.  These  calculations  are  based  onanozzle  section  one  inch 
long  located  at  the  throat.  Figure  22  depicts  the  rat''  at  which  heat  must  be 
removed  as  a  function  of  time  to  maintain  4800  F  wall  temperature.  The  dashed 
curve  in  the  same  figure  indicates  the  profile  achieved  by  varying  the  coolant 
flow  rate  linearly  starting  at  a  value  one-half  the  final  flow  rate.  Computer 
results  with  this  type  of  flow  rate  variation  are  presented  in  Section  2.3b  (3). 

(3)  Conduction  -  Convection  Computer  Study 

The  thermal  response  of  a  convect ively-cooied  heat  sink  was  studied 
in  more  detail  using  the  Lockheed  Thermal  Analyzer  program  described  in  Section  2.3a. 
The  analyses  to  date  have  been,  out  of  necessity,  a  study  of  t rends establ ished 
by  varying  heat  sink  thermal  conductivity,  heat  sink  size,  coolant  flow  rate  and 
composition,  hot  gas  side  convection  coefficient  and  nozzle  size  (scaling).  The 
exception  to  his  general  treatment  was  the  fairly  exact  analysis  of  the  sub-scale 
(1000  lbs.  thrust)  configuration  with  no  injection. 

Concurrent  with  this  computer  study,  the  determination  of  the  convection 
heat  transfer  coefficient  of  methane  including  the  effects  of  dissociation  has 
been  in  progress  (See  the  following  paragraph  4).  In  the  meantime,  the  values 
used  for  the  computer  input  were  calculated  using  the  properties  of  undissociated 
methane.  As  will  be  discussed  later,  the  actual  'alue  is  considerably  higher 
than  that  used  in  these  calculations.  This  will  mean  that  all  calculation  made 
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FIGURE  19.  COMPARISON  OF  TEMPERATURE  HISTORIES  WITH  AND  WITHOUT 
SURFACE  TEMPERATURE  CONTROL 
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FIGURE  21.  EFFECTS  OF  SURFACE  TEMPERATURE  CONTROL  ON  HEAT  TRANSFERRED  AND  STORED 
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FIGURE  22 .  COOLING  RATE  TO  MAINTAIN  CONSTANT  4800°F  WALL  TEMPERATURE 
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are  conservative.  Once  the  calculations  are  finalized,  a  previous  computer  run  will 
be  repeated  to  determine  the  effect.  Temperatures  measured  in  the  sub-scale  tests 
will  verify  these  effects. 

Figure  23  shows  typical  results  of  a  computer  run.  There  is  a  steep, 
rather  linear  temperature  gradient  between  the  surface  and  the  cooling  passage. 

From  the  cooling  passage  to  the  backvall  are  the  normal  transient  profiles.  For 
comparison,  a  temperature  profile  of  an  uncooled  heat  sink  of  the  same  configuration 
and  subjected  to  the  same  conditions  is  shown. 

A  search  through  the  literature  for  thermal  properties  of  pyrolytic 
graphite  indicated  that  the  only  widely  published  thermal  conductivity  data  at 
high  temperature  was  obtained  at  Raytheon  (Reference  10,  11,  12  &  13).  Another  source  of 
data  was  the  High  Temperature  Material  Inc.  catalogue  which  contained  extropol- 
ated  data.  The  latter  was  207.  lower  at  high  temperature  (4000  F  and  above). 

The  Raytheon  and  HTM  curves  are  shown  in  Figure  24  .  Each  of  these  curves 

were  input  into  the  computer  program  with  otherwise  indentical  conditions.  The 
following  table  shows  the  comparative  surface  temperatures  at  the  throat. 

TABLE  II 


Surface  Temperatures  for  Two  P.  G.  Conductivities 


Time(sec. ) 


Surface _ Temperature  (°F) 

HTM  Data  Raytheon  Data 


10 

4844 

4833 

20 

5165 

5094 

30 

5314 

5221 

40 

5404 

5305 

50 

5468 

5368 

60 

5515 

5418 

70 

5551 

5469 

80 

5580 

5493 

90 

5603 

5522 

100 

5622 

5547 

As  can  be  seen,  the  effect  on  wall  temperature  is  relatively  small. 
However,  as  the  surface  temperature  is  decreased  by  increasing  the  cooling  rates, 
the  difference  will  become  more  pronounced  because  of  the  higher  heat  flux. 

Since  the  rate  of  change  of  heat  storage  in  the  graphite  between  the  surface 
and  the  passages  is  small  (See  Figure  23),  a  steady  state  calculation  using  the 
equation 


q 


h  A 
g 


(T  -T  ) 
8  s 


k  A 
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can  be  made  to  determine  (T  -  T  )  for  any  values  of  k/x  and  T  .  Because  T 
cannot  be  determined  without  the^translent  solution,  only  the  relative  efflct 
of  k/x  can  be  determined.  It  should  be  noted  that  the  cooling  effectiveness  of 
convection  is  strongly  dependent  on  T  ,  hence  the  highest  possible  value  of  k/x 
is  deslreable.  Physically,  the  value^of  x  is  limited  to  approximately  0.25  or 
greater  for  designs  using  washers.  Because  pyrolytic  graphite  has  the  highest 
thermal  conductivity  of  any  suitable  materials  at  high  temperatures,  no  other 
material  is  currently  being  considered. 


A  computer  run  was  made  using  the  properties  of  heat  treated  pyroly¬ 
tic  graphite  obtained  from  Reference  13  .  The  conductivity  is  higher  than 

"as  deposited"  pyrographite  until  the  heat  treat  temperature  is  reached 
(approximately  5000  F).  The  wall  temperature  of  the  heat  treated  material 
lggged  that  of  the  "as  deposited"  material  somewhat  and  at  100  second  was  only 
9  F  lower.  Hence  the  results  indicated  that  heat  treatment  in  this  application 
was  not  warranted. 


Computer  runs  studying  the  effect  of  heat  sink  size  have  been  made, 
however,  optimization  studies  using  these  results  have  not  as  yet  been  made. 

The  study  will  involve  trade-offs  of  heat  sink  weight  vs.  coolant  system  weight. 

A  series  of  five  computer  runs  were  made  to  study  the  effects  of 
coolant  flow  rate  and  cooling  passage  diameter.  The  following  conditions  were 
fixed: 


T  -  6720  °F 

g 

h  *  0.00585  BTU/sec  in2  °F 

gt 

Dfc  -  2.50  in. 

D  (heat  sink)  *  6.00  in. 

X  -  0.25  in. 

Number  of  passages  *  50 

Coolant  -  methane 

k  -  Raytheon  Data  (Figure  24  ) 

The  variables  were  passage  diamete'  of  0.03  in.  and  0.04  in.  and  coolant  flow  rates 

of  0,  0.0325,  and  0.0475  lb/sec.  in  one  run  the  flow  was  varied  linearly  from 

0.0216  to  0.0434  lb/sec  with  an  average  of  0.0325  lb/sec.  The  effects  of  surface 
temperature  at  the  throat  is  shown  in  Figure  25 

Another  series  of  computer  runs  showed  the  effect  of  reducing  the  hot 
gas  side  heat  transfer  coefficient.  Since  the  actual  effects  of  film  injection 
is  as  yet  not  know,  the  runs  were  made  using  757*  and  50^  of  the  values  used  in 

other  runs  (See  Section  2.3c  ).  The  results  at  the  throat  are  shown  in  Figure  25  . 


-48- 


SURFACE  TEMPERATURE  (  F) 


1* 


6500 


6000 


5500 


5000 


4500 


4000 


UNCOOLED 


Wc-  0.0325  LB/ SEC 
D  -  0.040  IN. 


W  -  0.325  LB/SEC 
D  »  0.030  IN. 


.  W  -  0.0475  LB/SEC 
D  -  0.030  IN. 


in::: 


iiiiii  ******* 

:::::  n 


II  lilt!  HISS  ISttllS 

IMM11MMI  E*tt*«* 


!»:u»iuin;Eia»E; 

I«|I  hiimill  tilttMl' 


i:sri»ts;fEs;si  tnt 

Jiitriimtimii  tin  in 

!!!  !!!!«!!!!«!  (If  SHI1 


l||lll'r.|  ....  ""  rf 


ISHlit 


VARIABLE  W  WITH 
AVERAGE  •  0.0325  LB/ SEC' 
0  -  0.030  IN,  j 


0  20  40  60  80  100 


TIME  (SECONDS) 


R06702 
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Using  the  results  of  the  5000  lb.  thrust  nozzle,  a  rough  hand  calculation 
was  made  to  estimate  the  number  of  cooling  passages  and  the  methane  flow  rate  for 
tha  sub-scale  (1000  lb.  thrust)  nozzle  to  obtain  a  5000°F  surface  temperature  in 
100  sec.  The  result  indicated  28  holes  and  a  flow  rate  of  0.028  lb/sec. 

Runs  were  made  at  flow  rates  of  0.028  and  0.021  lb/sec.  The  surface 
temperatures  at  the  throat  at  100  second  wss  5280°P  and  5410  F  respectively.  The 
temperature  with  w  -  0.028  lb/sec  at  50  seconds  is  50?3°F.  Fifty  seconds  is  the 
planned  test  duration.  In  actuality,  this  temperature  will  probably  be  lower  since 
the  conservative  values  of  coolant  heat  transfer,  coefficient  were  used.  It  is 
interesting  to  note  that  the  5000  lb.  thrust  nozzle  Requires  0.16%  coolant  flow 
rate  compared  to  propellant  flow  rate  to  obtain  5300°F  surface  temperature, 
whereas  the  sub-scale  requires  0.50%  to  obtain  that  temperature. 

(4)  Methane  Convection  Coefficient 


As  shown  in  previous  analyses  for  100  second  firing  periods, 
the  coolant  rather  than  the  pyrolytic  graphite  mass  behaves  as  the  major  heat 
sink.  Therefore,  considerable  effort  must  be  devoted  to  formulation  of  analyti¬ 
cal  techniques  for  prediction  of  both  convective  and  radiative  heat  transfer  to 
the  coolant  from  adjacent  wall  surfaces.  As  previously  shown,  methane  gas  is  an 
attractive  coolant  from  materials  compatibility  and  protective  film  considerations. 
This  section  will  (1)  review  development  of  analytical  techniques  used  to  evaluate 
methane  as  a  heat  transfer  medium  and  (2)  present  results  of  a  typical  application 
of  methane  cooling  in  the  nozzle. 

The  laboratory  experiments  (Section  3.0  )  reveal  that  methane  disso¬ 

ciates  to  carbon  and  hydrogen  at  simulated  nozzle  cooling  passage  wall  temperatures. 
Carbon  particles  collected  from  these  tests  had  a  mean  diameter  of  about  0.1 
micron.  From  these  observations  it  is  apparent  that  a  meaningful  correlation 
for  convective  heat  transfer  to  the  methane  must  account  for  the  following: 

(1)  Solid  particle  production  and  growth  adjacent  to  the  heating 
surface. 

(2)  Methane  dissociation  near  the  heating  surface,  l.e.  effects  of 
chemical  reaction  and  composition  change. 

(3)  Selection  of  proper  reference  states  for  evaluating  pertinent 
properties  and  composition  of  the  coolant. 

Thermal  radiation  will  occur  between  the  wall  surface  and  coolant. 

The  following  steps  were  taken  to  evaluate  its  maximum  magnitude. 

(1)  Determination  of  gas  total  absorptivlties  for  black  body  radiation 
emitted  from  the  high  temperature  passage  surfaces. 

(2)  Determination  of  total  emissivlties  of  the  coolant. 

(3)  Calculation  of  maximum  radiating  heat  flux  and  equivalent  radiative 
heat  transfer  coefficient  (  h  )  for  comparison  with  convective 
heat  transfer  coefficients. 


-51- 


A  literature  survey  was  conducted  regarding  heat  transfer  to  fluids 
containing  suspensions  of  small  solid  particles.  It  is  consistently  reported 
that  presence  of  solid  particles  enhances  the  convective  film  coefficients. 

This  result  is  generally  interpreted  as  a  thickness  reduction  or  modif ication  of 
the  film  layer  adjacent  to  the  heating  surface  caused  by  the  turbulence  induced 
from  solid  particle  motion.  Most  information  relating  to  heat  transfer  to  gas- 
solid  mixtures  has  been  presented  since  1957.  At  this  time  Farber  and  Morley14 
presented  extensive  test  data  for  heat  transfer  to  an  alumina-silica  catalyst 
suspended  in  a  carbon  tetrachloride  fluid  system.  The  particle  size  (diameter) 
distribution  was  from  10  to  210  microns.  Farber  and  Morley  correlated  heat  trans¬ 
fer  data  for  the  mixture  flowing  in  a  tube  by  (Note:  Nomenclature  is  presented 
at  end  of  section): 


Nu  *  0.14  Re 


0.6 


0.45 


(1) 


As  a  typical  example  of  film  coefficient  increase  with  addition  of  solid  part¬ 
icles,  the  following  quantitative  observation  from  Far^gr  and  Morley's  data  is 
presented.  With  a  constant  gas  mass  velocity  of  1.50  ^  f  ®  the  convection 
film  coefficients  increased  by  factors  of  1.32  to  2.87  as  the  solids  loading 


ratios,  represented  by  Gg/G^. 


increased  in  value  from  2.0  to  10.5. 


The  correlation,  equa* ir  1,  was  derived  for  specific  fluids,  solids 
and  particle  size  distributions,  ar.u  therefore  is  probably  not  suitable  for  the 
dissociating  methane  system.  In  fact,  it  appears  that  particle  size  is  a  8overn* 
ing  parameter  in  the  convective  heat  transfer  process.  Wen  and  Miller  very 
roughly  correlated  the  convective  heat  transfer  for  suspensions  of  solids  with  an 
equation  containing  particle  diameters  as  a  variable.  It  was  also  very  recently 
suggested^  that  particle  size  and  composition  would  be  of  importance  in  a 
general  type  correlation.  There  is,  however,  no  acceptable  general  correlation 
in  existence  at  present. 


Very  fortunately  there  is  heat  transfer  data  available  for  solid 
suspensions  with  particle  sizes  and  compositions  similar  to  those  anticipated 
in  the  dissociating  methane.  In  1959  and  1960  the  Babcock  and  Wilcox  Company 
presented  results  of  an  extensive  heat  transfer  and  pressure  drop  testing 
program  conducted  with  graphite  suspensions  in  nitrogen,  helium,  and  Freon 
(CF  )  gas  flows.  The  primary  goal  of  the  program  was  to  develope  means  for 
increased  heat  transfer  performance  in  gas  cooled  nuclear  reactors.  In 
general,  the  graphite  particle  diameters,  as  determined  by  electron  micro¬ 
scope  examination,  varied  between  0.1  and  6  microns.  Mean  particle  diameters 
were  approximately  1  micron.  The  correlation  established  from  the  Babcock  and 
Wilcox  heat  transfer  tests  in  smooth  tubes  is: 
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As  shown  in  Reference  19, this  correlation  is  applicable  for  a  wide  range  of 
graphite  loadings  of  nitrogen,  helium  and  CF,  gases.  The  subscript  g  signifies 
that  the  parameter  (RePr)  Is  evaluated  for  tne  gaseous  component  of  the  flow 
only.  It  Is  assumed  that  pertinent  properties  are  evaluated  at  bulk  gas  temper¬ 
atures  . 

In  the  present  case  of  cooling,  the  methane  operates  over  a  greater  wall 
to  bulk  fluid  temperature  difference.  It  was  therefore  judged  necessary  to  re¬ 
establish  a  correlation  Interms  of  parameters  to  be  evaluated  at  specified 
reference  temperatures.  The  correlation  is  developed  from  the  basic  Babcok  and 
Wilcox  heat  transfer  data.  Presented  In  Reference  19  ,  this  data  is  specifi¬ 

cally  for  flow  consisting  of  various  solids  loadings  of  graphite  In  nitrogen, 
helium,  and  CF^.  Solid  to  gas  weight  fractions  ranged  from  1.14  to  10.33  In  the 
nitrogen,  8.65  to  81.4  In  helium  and  0.360  In  CF^.  All  tests  were  conducted  In  a 
0.313  Inch  I.D.  tube.  Oelssler  and  Presler  ®  conducted  analysis  on  the  determin¬ 
ation  of  reference  temperatures  for  property  evaluation  for  heat  transfer  to 
helium  ,  argon,  hydrogen,  air,  and  carbon  dioxide  under  conditions  of  wall  temp¬ 
erature  variation  from  2000  to  6000  F.  The  recommended  reference  temperature  Is 
calculated  from  T  »  T  +  x  (T  -  T^).  The  quantity  x  is  a  function  of  the  fluid 
Reynold's  Number  2nd  tne  wall  temperature.  Specific  reference  Is  made  to  Figure  27. 
• 

A  correlation  of  the  form 


Nu  -  C  Re  °*8Pr  0,4 

gx  1  gx  gx 


W  C  A  . 

S  8  +  1 

w  c 

g  8 
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was  established  from  the  Babcock  and  Wilcox  data  by  evaluating  gas  properties 
at  appropriate  reference  temperatures.  Figure  2b  illustrates  the  correlation 
of  various  data  points  and  establishes  the  constant  C  *  0.023.  As  In  the 
Babcock  and  Wilcox  correlation,  Equation  2,  the  Nusseit,  Reynolds  and  Prandtl 
numbers  are  evaluated  from  the  gas  properties  only. 

In  the  case  of  dissociating  methane  allowance  must  be  made  for  effects 
of  chemical  reaction  on  heat  transfer  rates.  As  explained  by  Butler  and  Brokaw?l, 
reacting  gases  In  chemical  equilibrium  may  have  considerably  higher  effective 
conductivities  than  gases  of  frozen  composition.  The  Increase  In  the  conductivity 
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FIGURE  27.  VALUES  OF  X  FOR  CALCULATING  REFERENCE  TEMPERATURES  FROM 
DEISSLER  AND  PRESLER 
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is  interpreted  as  a  result  of  transport  of  chemical  enthalpy  of  molecules  which 
diffuse  according  to  concentration  gradients  established  by  chemical  reaction.  The 
concentration  gradients  and  therefore  the  effective  thermal  conductivities  are 
governed  by  the  variation  in  the  equilibrium  gas  composition  with  temperature. 

Considerable  information  is  available  on  suitable  convective  heat  trans¬ 
fer  correlations  for  flow  with  chemical  reactions.  Relatively  low  velocity  turbu¬ 
lent  pipe  flow  heat  transfer  with  chemical  reaction  has  been  correlated  by  the 
following  approach  according  to  Bartz  1  : 


q  * 


S'(^u  )  (hw 


The  Stanton  number  is  identical  to  that  defined  for  non-reactlve  flow.  As 
indicated,  the  enthalpy  difference  of  the  fluid  at  wall  and  bulk  states  is 
essentially  the  driving  potential  for  heat  transfer.  Additional  background  for 
the  general  use  of  equation  4  for  chemically  reactive  heat  transfer  is  presented 
in  References  25  and  26  by  E.  R.  G.  Eckert.  The  major  restriction  on  its 
use  is  the  requirement  of  unity  Lewis  Number.  Aa  stated  in  Ibele27  and  interpret¬ 
ed  from  Altman  and  Wi„e  28  ,  the  Lewis  number  very  closely  approaches  unity  for 
turbulent  flow. 


Turbulent  Lewis  number  ■  t  ■  1  because  t-  ;>P, 

and  magnitudes  of  D  and  oC are 
not  too  different 

where  D  ■  moleculardif fusivity  ^ 

o<  -  thermal  dlffuslvlty,  “7^ — 

C-  »  turbulent  eddy  d  if  fusivity  ^p 

The  methane  flowing  through  the  nozzle  cooling  passages  will  exhibit  considerable 
turbulence.  Based  on  gas  phase  properties  and  tube  diameter  the  Reynold's  Numbers 
are  within  the  range  5,000  to  11,000.  Furthermore,  the  formation  and  growth  of 
carbon  particles  will  induce  turbulence  near  the  wall-fluid  interface. 

The  correlation,  equation  3  for  solid  suspension  flow  is  simply  changed 
into  the  form  of  equation  4: 
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from  Equations  3  and  5: 
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from  Equations  4  and  6: 
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It  is  again  pointed  out  the  subscript  g  refers  to  parameter  and/or 
property  evaluation  for  the  gaseous  components  only. 


In  application  of  Equation  7.  the  selection  of  a  proper  reference  state 
for  property  ev  luation  poses  some  question.  If  the  assumptions  leading  to  the 
reference  temperature  technique  presented  in  Reference20  remain  valid  in  the  case 
of  dissociating  flow  with  solid  particle  production,  then  Equation  7  properties 
may  be  evaluated  accordingly.  An  alternative  method,  which  may  prove  applicable, 
has  been  presented  by  Eckert^-*  ,  in  which  properties  are  evaluated  at  a  reference 
enthalpy.  This  approach,  however,  is  especially  adapted  for  heat  transfer  from 
high  velocity  flows. 


As  an  example  of  the  application  of  the  derived  correlation,  Equation  7, 
heat  flux  and  corresponding  convective  film  coefficients  were  calculated  for  com¬ 
binations  of  wall  temperature  and  bulk  coolant  temperature  expected  to  exist  at 
four  positions  along  the  nozzle  cooling  passages.  Results  are  presented  on  Figure 
29  .  For  purposes  of  comparison,  coefficients  calculated  for  the  case  of  non¬ 

dissociating  methane  gas  are  also  included.  Physical  properties  in  both  cases 
were  evaluated  at  reference  temperatures  according  to  Reference  20  and  shown 
in  Figure  27  .  For  reference  temperature  selection,  coolant  composition 

adjacent  to  the  heating  surface  was  assumed  to  be  primarily  hydrogen.  In  the 
comparison  case  of  non-dissociating  methane  gas,  the  Reynold's  number  was  based 
on  bulk  fluid  mass  velocity.  The  pertinent  properties  and  parameters  required 
in  these  calculations  are  shown  in  Figures  30,  31,  and  32.  The  specific  heat, 
thermal  conductivity,  and  viscosity  data  for  methane  and  hydrogen  were  obtained 
from  Svehla^ .  Properties  of  equilibrium  hydrogen  and  methane  mixtures  were 
then  determined  from  alignment  charts  presented  in  Brokaw^O .  Equilibrium  com¬ 
position  of  methane  and  its  dissociating  products  were  obtained  from  an 
Aeronutronic  computer  program  in  which  JANAF  Thermochemical  Data  ^1  t8  utilized. 
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FIGURE  29.  COMPARISON  OF  ESTIMATED  CONVECTIVE  FILM  COEFFICIENTS 
FOR  DISSOCIATING  AND  NON -DISSOCIATING  METHANE 
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FIGURE  30.  VISCOSITY,  THERMAL  CONDUCTIVITY,  AND  SPECIFIC  HEAT 
OF  METHANE  AND  GASEOUS  DISSOCIATION  PRODUCTS 


Xi/aD  SVD  iO 


f  WgV, 

\WgCg 


FIGURE  31.  DENSITY,  MASS  FRACTION,  AND  FACTOR 


OF  METHANE  AND  GASEOUS  DISSOCIATION  PRODUCTS 
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ENTHALPY  OF  METHAlt  AND  ITS  GASEOUS  DISSOCIATIO 
PRODUCTS  AT  700  PSIA 


The  enthalpy  term  of  Equation  7  refers  to  enthalpies  of  gaseous  components  only 
by  virtue  of  the  definitions  of  demensionless  parameters  of  equations  3  and  5. 

These  enthalpies,  shown  on  Figure  32  ,  are  also  obtained  from  the  Aeronutronic 
computer  program  and  JANAF  Thermochemical  Data. 

The  information  presented  on  Figure  29  shows  that  the  convective  film 
coefficients  with  dissociation  and  solid  carbon  suspension  is  roughly  twice  that 
calculated  for  heat  transfer  to  non>reacting  methane  for  identical  passage  wall  and 
bulk  fluid  temperatures  levels. 

In  an  effort  to  obtain  the  magnitude  of  thermal  radiation  an  analysis 
has  been  made  at  the  location  of  highest  passage  wall  temperature.  Radiation 
properties  of  the  coolant  were  obtained  from  a  paper  by  Stull  and  Plass  ^  which 
presents  monochromatic  emissivity  data  for  thermally  decomposing  hydrocarbon 
gases  containing  various  sizes  and  distributions  of  carbon  particles.  According 
to  this  data,  the  decomposing  methane  in  the  highest  temperature  region  of  the 
nozzle  cooling  passages  will  behave  as  a  black  body  over  the  entire  range  of 
applicable  radiation  wavelengths.  With  the  assumption  that  the  pyrolytic 
graphite  cooling  passage  wall  behaves  as  a  black  body,  a  net  radiative  heat 
flux  may  be  determined  by  Equation  8: 


q  - 


(8) 


An  equivalent  radiative  heat  transfer  coefficient  may  be  found  from 
Equation  9. 
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o 

For  the  maximum  wal^  surface  temperature  of  3470  F  and  corresponding  bulk  coolant  ^ 
temperature  of  1414  F  the  radiative  heat  transfer  coefficient  is  equal  to  3.67x10 

_ BTJL _  _ This  is  far  less  than  the  convective  heat  transfer  coefficient  which 

sec  in  oR. 

has  a  magnitude  of  78.4x10  BTU _  at  the  same  location. 

.  2  °B 
sec  in  R 

The  following  observations  summarize  the  major  results  of  a  heat 
transfer  analysis  to  the  nozzle  cooling  passages: 

1.  Dissociating  methane  coolant  yields  higher  convective  film  heat 
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transfer  coefficients  than  non-dissociating  methane.  The  increase 
is  due  primarily  to  the  chemical  reaction  itself.  The  effect  of  small 
solid  carbon  particle  formation  in  increasing  convection  is  secondary 
according  to  correlations  developed  from  data  of  References  17,  18  61  19. 

2.  Radiation  heat  transfer  from  the  wall  surface  to  the  coolant  is 

lower  than  convective  heat  transfer,  though  it  was  determined  that 
both  coolant  and  wall  surface  behave  as  black  bodies. 

Additional  work  will  be  required  in  the  following  areas  in  order  to  complete  the 
analytical  treatment  of  cooling  passage  heat  transfer. 

1.  Thorough  examination  of  applicability  of  various  reference  temperatures 
or  enthalpies  for  property  evaluation  in  the  correlation  for  convective 
heat  flux.  The  importance  of  reference  point  selection  is  illustrated 

in  Figures  30  and  31  where  properties  vary  considerably  with  temperature. 
The  Deissler  and  Presler^Oreference  temperature  technique,  which  is 
used  to  predict  results  of  Figure  29,  is  derived  originally  for  non- 
reacting  flows.  Should  an  alternative  approach  be  required,  a  modifi¬ 
cation  of  the  Eckert  reference  enthalpy  parameter  may  be  utilized. 

Heat  transfer  data  from  the  Initial  subscale  tests  will  be  useful  in 
evaluating  selection  of  reference  points. 

2.  The  objective  of  the  nozzle  coolant  heat  transfer  analysis  is  to 
obtain  input  data  for  the  thermal  analyzer  digital  computer  program 
which  predicts  temperature  profiles  through  the  pyrolytic  graphite 
heat  sink  and  coolant  bulk  temperatures.  Because  the  convective 
heat  transfer  coefficients  are  functions  of  both  wail  and  bulk 
coolant  temperatures,  it  will  be  necessary  to  establish  an  Iteration 
procedure  in  the  computer  program. 
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NOMENCLATURE 


c 

P 

d 

G 

h 

H 
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Nu 

Pr 

q 

Re 

St 

T 

U 

V 

/' 


o  •  • 

Subscripts 

b  .  .  . 

g  •  •  • 

at  . 
r  . 

w  . 

x  .  . 


BTU 

Specific  heat, 


Tube  diameter,  in. 


Mass  velocity, 


lb 


sec  in 


Convective  film  coefficient 


BTU 


.  2  o 

sec  in  F 


Enthalpy,  BTU/lb 
Thermal  conductivity, 


Nusselt  number,  ~ 

k 


BTU 


,  2  °* 
sec  in  F 


in 


Prandtl  Number, 


/U 


Heat  flux,  BTU/sec  in 


Reynold's  Number, 


pUd 


Z4 


Stanton  Number 


U  Cp 


Temperature,  °F,  (°R  for  radiation) 
Velocity,  in/sec 
Flow  rate,  lb/sec 
lb 

Viscosity, 

Density , 

Stefan-Bol tzmann  Constant 


lb 

in- 


Bulk  fluid  states  or  properties 
Denotes  gaseous  components 
Denotes  mixture 

Pertaining  to  thermal  radiation 
Wall  states  or  properties 
Reference  point  states  or  properties 
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Flame  Side  Convection 


c . 


The  primary  objective  of  the  film  injection  is  to  modify  the  chemical 
composition  of  the  boundary  layer  in  the  nozzle  throat  region  in  such  a  way  that 
chemical  attack  is  reduced  or  eliminated.  Because  the  convection  fluid  cannot  be 
injected  at  temperatures  which  approach  the  design  surface  temperature,  a  signi¬ 
ficant  film  cooling  effect  and  consequently  an  increase  in  coolant  thermal  effi¬ 
ciency  will  result.  The  investigation  of  both  the  chemical  and  heat  transfer 
effects  requires  a  complete  characterization  of  the  convection  boundary  layer 
without  injection.  Such  analysis  also  provides  the  basis  for  the  design  of  the 
convectively  cooled  (no  injection)  heat  sink  nozzle.  It  should  also  be  recogni¬ 
zed  that  the  nozzle  convective  heat  transfer  problem  remains  the  subject  of  general 
interest  within  the  industry.  During  the  first  quarter,  efforts  were  devoted  to 
(1)  surveying  the  recent  literature  pertaining  to  boundary  layer  analysis,  (2) 
evaluating  the  J.P.L  nozzle  heat  transfer  program,  (3)  performing  preliminary  heat 
transfer  calculations  for  the  sub-scale  nozzle  designs,  and  (4)  the  selection  of 
an  improved  radiation  heat  transfer  analysis. 

Recently  reported  results  of  nozzle  heat  transfer  studies  have  provided 
some  additional  enlightenment  but  no  major  progress  has  been  recorded.  The  exist¬ 
ence  of  laminar  and  transitional  boundary  layer  flows  in  nozzles  has  been  demon¬ 
strated  in  Reference  33  .  From  these  results  it  is  evident  that  turbulent  bound¬ 
ary  layers  will  persist  throughout  the  inlet  and  throat  sections  of  nozzles  to  be 
tested  in  this  program.  This,  of  course,  is  not  an  unexpected  result  for  gas 
flows  and  the  presence  of  a  condensed  phase  is  expected  to  increase  turbulence. 
Reference  33  has  also  made  contributions  to  the  understanding  of  the  effects  of 
boundary  layer  acceleration  and  nozzle  entrance  conditions.  It  has  also  become 
apparent,  from  the  literature  discussed  in  Section  2.2  that  the  effects  of  trans- 
pirational  injection  of  resin  pyrolysis  products  in  the  chamber  and  upper  inlet 
section  may  have  a  significant  influence  on  boundary  layer  development.  However, 
it  seems  likely  that  such  effects  could  easily  be  obscured  by  acceleration, 
radiation,  and  condensed  phase  effects.  It  can  also  be  deduced  from  Reference  33 
that  the  anticipated  variations  of  chamber  and  nozzle  surface  temperatures  will 
have  a  significant  influence  on  boundary  layer  development  in  addition  to  the 
usual  effects  on  radiation  and  gas  properties. 

Considering  all  the  aspects  of  the  boundary  layer  development  and  heat 
transfer  problem,  the  most  logical  approach  to  predicting  heat  transfer  is  to  use 
existing  methods  in  such  a  way  that  the  results  can  confidently  be  expected  to  be 
conservative.  Then,  an  overcooled  nozzle  design,  when  properly  Instrumented,  can 
be  expected  to  provide  reasonable  empirical  corrections  to  the  original  predictions. 
Data  obtained  from  undercooled  nozzle  tests  is  usually  compromised  by  partial  or 
total  failure  of  the  inlet  or  throat  surfaces.  The  convection  cooled  nozzle  de¬ 
signs  (no  injection)  are  expected  to  provide  the  means  to  make  empirical  estimates 
or  corrections  to  the  predicted  heat  transfer.  It  should  be  noted  that  the  term 
overcooled  really  implies  that  superior  structural  and  erosion  performance  is  the 
observed  result  of  actual  testing.  This  is  simply  because  the  maximum  or  optimum 
design  operating  conditions  cannot  as  yet  be  predicted  analytically. 

At  the  start  of  this  program  it  was  decided  to  use  the  modified  Barts 
method  for  calculating  boundary  layers.  The  method  and  associated  computer  pro¬ 
gram  are  well  described  in  Reference  34  and  will  not  be  presented  here.  This 
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analysis  has  been  reviewed  to  determine  its  applicability  for  boundary  layer  cal¬ 
culations,  both  with  and  without  injection.  It  was  found  that  a  number  of  minor 
improvements  could  be  made  to  both  the  theory  and  programming.  These  have  not 
been  carried  out,  primarily  because  the  program  may  also  ultimately  be  revised  for 
use  with  the  film  analysis.  For  the  no  injection  case,  the  major  difficulty  in 
applying  the  analysis  arises  from  the  spatial  and  temporal  variations  in  surface 
temperature  (or  enthalpy  of  the  gas  at  the  wall).  Fortunately,  the  convection 
cooled  nozzle  design  operates  with  relatively  constant,  in  time,  surface  tempera¬ 
tures  and  the  important  axial  variations  in  surface  temperature  can  be  predicted 
without  extensive  iteration,  however,  it  is  important  to  realize  that  the  heat 
transfer  during  the  major  surface  temperature  transient  (the  first  20  to  30 
seconds  of  operation)  cannot  be  easily  predicted.  Such  calculations  would  be 
possible  if  the  behavior  of  the  alumina  were  known. 

For  the  film  injection  case,  the  computer  program  could  be  used,  to 
estimate  film  heat  transfer  and  chemical  mixing.  This  method  was  not  described 
in  Section  2.2  and  is  referred  to  as  a  far-field  solution.  It  would  be  assumed 
that  the  momentum  mixing  region  is  negligible.  Thus,  a  new  1 /7th  power  profile 
could  be  constructed  at  the  point  of  injection  via  the  conservation  equations. 

The  initial  wall  temperature  would  be  taken  to  be  equal  to  the  fluid  temperature 
at  the  point  of  injection.  A  simplified  diffusion  analysis  could  then  be  formu¬ 
lated  to  determine  the  chemical  mixing  and  the  consequent  density  variations.  Pre¬ 
liminary  considerations  suggest  that  this  technique  may  have  merit  if  more  accurate 
solutions  appear  to  be  impractical. 

Several  computer  calculations  were  performed  to  provide  a  preliminary 
estimate  of  the  nozzle  heat  trasnfer.  The  predictions  were  slightly  higher  than 
the  values  deduced  during  contract  AF  04(611 ) -8387  for  the  same  propellant.  To 
facilitate  the  calculations,  on  arbitrary  geometry  was  chosen  and  surface  tempera¬ 
tures  were  held  constant  at  4300  F.  The  results  of  these  calculations  are  not 
presented  in  detail  because  neither  the  geometry  nor  temperature  conditions  have 
been  preserved  in  the  final  designs.  Improved  geometry,  property  data,  enthalpy 
and  surface  temperature  data  will  be  utilized  in  subsequent  computer  calculations. 

During  the  next  quarter  of  the  contract  period,  a  radiation  analysis 
of  the  convergent  section  will  be  performed.  The  stimulus  for  this  analysis  was 
that  recent  data  obtained  at  Philco  Research  Laboratories  on  the  emissivity  of  a 
two  phase  mixture  indicates  that  the  radiation  heat  flux  may  be  somewhat  higher 
than  previously  suspected.  Presented  below  are  the  general  expressions  that  will 
be  used  '  r  analysis  of  the  gaseous  -par  tic  le  radiation  heat  flux.  Using  the  network 
method,  the  following  equations  can  be  written: 

VAi  -  <Ei  •  V 

Vi  ■  ViAi +  <>  -  6i>¥ 
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where : 


=  heat  flux  to  A^ 

-  emissivity  of  area  A^  under  consideration 
E^  =  emissive  power  of 
*  radiosity  of  A^ 

G  *  irradiation  to  A^ 


Assuming  that  the  particles  may  be  treated  as  a 
the  irradiation  from  the  two  phase  mixture  and  an  area  A 


35 


and 


36 


gas  (including  scattering), 
to  an  area  A^  is  References 


AjG 


J2 


cos  /3 ^  cos  /3 , 

^  2 

77  r 


dA2  dAl 


+ 


cos^Jg  d/3g 


dA 


1 


where : 


E  *  emissive  power  of  two  phase  mixture 


transmissivity  and  emissivity  of  two  ph“se  mixture 

angle  between  normal  of  A  and  line  connecting  centers 
of  Ax  and  A2  1 

angle  between  normal  of  A2  and  line  connecting  centers 
of  A^  and  A2 


azlmth  angle 
polar  angle 


It  may  be  noted  that,  for  two  phase  gaseous  radiation,  one  employs 
spherical  coordinates,  as  the  transport  of  energy  is  a  volume  phenomenon.  With 
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particles,  the  transmissivity  and  emissivity  of  a  gas  are  defined  as 


6  -  1  -  Z  -  1  -  e_as 

where : 

s  a  path  length  ' 

a  -  absorption  coefficient  including  scattering  induced 
by  the  particles. 


At  Philco  Research  Laboratories,  experiments  are  being  performed 
[under  contract  NOnr  3907(00)]  on  emissivity  measurements  of  an  A1  0.  -HO 

L  j  2 


rocket  engine.  The 
to  those  that  exist 
will  be  utilized  in 


particle  concentration,  size  and  temperature  are  equivalent 
in  the  film  protected  nozzle.  The  results  from  these  experiments 
the  radiation  analysis. 
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SECTION 


3 


LABORATORY  EXPERIMENTS 


3 . 1  CHEMICAL  REACTIONS 

A  series  of  laboratory  experiments  have  been  initiated  to  define  the 
extent  of  interaction  of  the  various  promising  injection  gases  with  the  graphite 
capillaries.  An  attempt  has  been  made  to  simulate  the  nozzle  design  conditions 
of  flow  velocity,  temperature  gi*dient  and  pressure.  Three  grades  of  graphite 
were  tested,  Great  Lakes  HLM-85,  Basic  Carbon  Graph-i-tlte  GX  and  pyrolytic  graph* 
ite. 


a.  Experimental 


The  experimental  apparatus  is  shown  schematically  in  Figure  33 
Essentially,  it  consists  of  all  stainless  steel  high  pressure  chamber  that  can  be 
evacuated  and  then  back-filled.  The  system  has  been  designed  to  accommodate  gas 
pressures  to  600  pslg.  All  tests  were  performed  under  forced  gas  flow  conditions 
with  provision  for  metering  of  gas,  and  ballast  valves  for  maintaining  constant 
pressure. 


Test  specimens  consist  of  graphite  rods,  7/8-lnch  diameter  by  8-1/2 
Inches  long,  with  the  top  2-1/2  Inches  containing  the  capillaries  to  be  tasted. 

The  specimens  are  inductively  heated  by  means  of  a  Lepel  vacuum  tube  radio  fre¬ 
quency  generator  having  an  operating  frequency  of  450  kilocycles  and  an  output 
rating  of  20KW.  The  furnace  consists  of  a  multi-turn,  barrel-shaped  coll  wound  in 
a  helix  around  the  test  specimen. 

Operational  procedure  Involves  first  pressurising  the  chambmr  to  250 
psla  with  argon  and  then  Introducing  the  injection  gaa  through  the  capillary  at 
a  predetermined  flow  rate  (STP)  and  a  constant  delivery  pressure  of  350  pale 
(differential  pressure  cf  100  psla).  After  passing  through  the  capillary,  the 
Injection  gas  mixes  freely  with  the  argon  in  the  chamber.  Constant  chamber 
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pressure  is  maintained  by  the  use  of  micrometer  ballast  valves  at  the  chamber  and 
at  the  gas  supply  cylinder. 


Sample  temperatures  are  measured  with  a  Pyrometer  Instrument  Company 
disappearing  filament-type,  optical  pyrometer  by  sighting  on  the  exterior  surface 
of  the  specimen  through  a  ground  quartz  window.  Thus,  all  recorded  temperatures 
are  brightness  temperatures  of  the  exterior  surfaces.  It  is  recognized  that  the 
capillary  wall  temperatures  may  be  one  hundred  to  two  hundred  degrees  cooler. 

All  tests  were  performed  at  an  exterior  surface  tempe.ature  of  2100  C. 

The  power  is  turned  on  immediately  after  injection  gas  flow  is  started. 
Approximately  15  seconds  are  required  to  bring  the  graphite  sample  to  the  test 
temperature  of  2100  C.  Where  possible,  tests  ran  for  three  minutes.  Shorter 
duration  tests  were  necessary  in  some  of  the  methane  pyrolysis  experiments  due  to 
accumulation  of  carbon  on  the  furnace  coil  with  resultant  electrical  arcing. 

b.  esults 


(1)  HLM-85 


Six  gases;  methane,  carbon  dioxide,  ammonia,  boron  trifluoride, 
carbon  monoxide,  and  hydrogen,  were  tested  for  their  compatibility  with  HLM-85 
graphite  at  or  near  2100  C  by  passing  the  injected  gas  through  a  0.040-inch  dia¬ 
meter  by  2.5  inches  long  capillary.  With  the  exception  of  ammonia,  gas  flow 
rates  varied  from  60  to  120  cubic  feet  per  hour.  Ammonia  is  available  commerci¬ 
ally  only  as  a  liquified  gas  under  its  own  vapor  pressure  of  130  psia  at  70  F. 

In  order  to  maintain  a  differential  pressure  of  100  psia,  the  chamber  was  pressur¬ 
ized  only  to  30  psia  with  argon.  Under  these  conditions,  ihe  maximum  flow  rate 
which  could  be  achieved  with  ammonia  was  8  cubic  feet  per  hour. 

The  results  are  shown  in  Figures  34  through  37  .  All  tests  were  of 

three  minuces  duration.  Figure  34  represents  a  side  view  of  the  test  speci¬ 

mens  after  test  and  Figure  35  a  top  view.  Major  chemical  attack  of  the  graphite 
capillary  was  noted  for  ammonia  and  CO^  with  probable  formation  of  gaseous  CN 
and  CO,  respectively,  as  products  of  reaction.  Minor  corrosion  (  10%  hole  en¬ 

largement)  was  noted  for  BF^  and  H^.  In  the  case  of  methane,  large  quantities  of 
amorphous  carb  n  were  found,  not  only  on  the  exterior  surface  c£  the  test  speci- 
ment  as  seen  in  Figure  34  ,  but  also  over  the  entire  chamber  assembly.  The 
0.040-inch  capillary  was  found  to  be  completely  open,  with  no  evidence  of  carbon 
deposition  along  the  wails.  Figures  36  and  37  are  cross-sections  of  the 
specimens  after  test.  These  were  prepared  by  carefully  machine  milling  through 
the  specimens  until  the  capillary  walls  were  exposed.  The  material  around  the 
CH^  test  specimen  is  an  organic  "potting"  resin,  used  to  insure  that  the  amorphous 
carbon  on  the  exterior  of  the  specimen  did  not  flake  off. 

Figure  38  is  an  electron  photomicrograph  (10000  X)  of  the  carbon 
products  collected  from  the  top  area  near  the  orifice  of  the  CH^  test  specimen. 

The  carbon  particles  were  spherical  and  ranged  in  size  from  0.1  to  0.05  microns. 

As  can  be  seen,  the  particles  were  highly  agglomerated. 
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FIGURE  34.  HLM-85  GRAPHITE  -  APPROXIMATELY  2/3X 


R06522 


HLM-85  GRAPHITE  -  TOP  VIEW  -  APPROXIMATELY  2/3X 
ORIGINAL  HOLE  SIZE  -  0.040  INCH 


FIGURE  35. 


R06S20 


FIGURE  37.  HLM-85  GRAPHITE  -  SECTION  VIEW  AFTER  TEST  - 
APPROXIMATELY  1.2X 


FIGURE  38.  ELECTRON  PHOTOMICROGRAPH  OF  CARBON  PARTICLES  IOOOOX 
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(2 )  Graph- 1 -tite 


Since  no  evidence  of  carbon  deposition  along  the  capillary  walls 
was  found  when  methane  was  injected  over  1ILM-85,  the  test  was  repeated  with  Basic 
Carbon  Graph-i-tite  GX,  a  higher  density,  finer  grained  graphite.  The  results  for 
a  three  minute  test  at  a  methane  flow  rate  of  60  cubic  feet  per  hour  were  essenti¬ 
ally  the  same  as  for  HLM-85.  No  evidence  of  carbon  deposition  along  the  capillary 
walls  was  found.  Considerably  less  carbon  was  found,  however,  along  the  exterior 
walls  of  the  test  specimen,  probably  reflecting  the  decrease  in  porosity  of  Graph- 
i-tite  over  HLM-85.  Carbon  particles  collected  from  the  top  area  near  the  ori¬ 
fice  were  spherical  and  ranged  in  size  from  0.1  to  0.05  microns.  Again,  a  high 
degree  of  agglomeration  was  noted. 

Since  the  carbon  particles  presumably  originate  in  the  capillaries  from 
the  pyrolysis  of  methane,  the  question  arose  as  to  whe  .her  or  not  sharp  angles  in 
the  Graph-i-tite  capillary  walls  would  cause  deposition,  i.e.,  will  the  carbon 
particles  follow  the  main  gas  stream  in  making  sharp  bends  or  will  they  accumulate 
at  the  bend.  In  order  to  study  this,  the  capillary  was  modified  to  include  two 
135  degree  turns,  as  shown  in  cross-section  Figure  39  .  The  entrance  capillary, 

A,  was  0.100-inch  diameter  by  2  inches  .ong;  and  the  two  exit  capillaries,  B  and 
C  were  0.040-inch  diameter  by  9/16  inch  long.  The  junction,  D,  thus  represents  a 
plenum  where  carbon  particle  accumulation  would  be  expected.  Two  tests  were  run, 
one  with  a  total  methane  flow  rate  of  approximately  120  cubic  feet  per  hour  (60 
cubic  feet  per  hour  per  exit  capillary)  and  the  second  with  a  total  methane  flow 
rate  of  approximately  60  cubic  feet  per  hour  (30  cubic  feet  per  hour  per  exit 
capillary).  The  first  speciment  (120  cubic  feet  per  hour  methane  flow)  was  run 
at  temperature  for  75  seconds  and  the  second  for  60  seconds.  In  both  cases,  the 
experiments  were  terminated  due  to  accumulation  of  carbon  on  the  furnace  coil  with 
resultant  electric  arcing.  As  can  be  seen  in  Figure  39  ,  no  accumulation  of 

carbon  was  found  in  either  test.  Figures  40  and  41  are  side  views  of  the  first 
specimen  (120  cubic  feet  per  hour  methane  flow)  showing  the  two  open  capillaries. 

It  is  of  interest  to  note  that  the  grey  pyrolytic  graphite  areas  appear  to  reflect 
the  exit  methane  flow  patterns.  From  electron  microscope  studies,  carbon  particles 
near  the  exit  holes  were  again  found  to  be  spherical  and  ranged  in  size  from  0.1 
to  0.05  microns. 

An  additional  experiment  was  made  with  Graph-i-tite  in  which  two 
"straight-through"  capillaries,  each  0.040-inch  diameter  by  2-1/2  Inches  long, 
were  provided.  The  total  methane  flow  rate  was  approximately  120  cubic  feet  per 
hour  (60  cubic  feet  per  hour  per  exit  capillary).  Again,  no  accumulation  of  carbon 
was  found  in  either  capillary  as  evidenced  in  Figure  38. 

(3)  Pyrolytic  Graphite 

Pyrolytic  graphite  washers  were  stacked  in  an  arrangement  shown 
schematically  in  Figure  42  .  The  ATJ  graphite  shell  served  as  a  suscept-heater . 
Two  tests  were  run  ,  both  with  a  methane  flow  of  approximately  60  cubic  feet  per 
hour.  In  th<»  first  test,  lasting  70  seconds,  the  top  washer  shifted  slightly 
during  the  run,  effectively  closing  the  exit  capillary.  As  a  result,  the  carbon 
particles  from  the  pyrolysis  of  the  methane  backed  up  and  clogged  the  capillaries 
of  the  remaining  pyrolytic  graphite  washers.  For  the  second  test,  the  top  washer 
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FIGURE  39.  GX  GRAPH- 1 -THE  -  AFTER  TEST  APPROXIMATELY  1 . 3X 


FIGURE  40.  GX  GRAPH-I-TITE  -  AFTER  TEST  APPROXIMATELY  1.5X 
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FIGURE  41.  GX  OKAPH-I-TITE  -  AFTER  TEST  APPROXIMATELY  1.5X 
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FIGURE  42. 


PYROLYTIC  GRAPHITE  WASHER  SPECIMQ*  -  SCHB1ATIC 


was  securely  anchored  to  the  ATJ  graphite  sleeve  by  three  dowel  pins.  This  re¬ 
sulted  in  a  satisfactory  test  lasting  for  60  seconds.  After  test,  the  capillary 
was  found  to  be  completely  open  with  no  evidence  of  carbon  deposition  along  the 
capillary  walls. 

c .  Conclusions 

(1)  Ammonia  and  carbon  dioxide  severely  attack  graphite  at  the  tempera¬ 
tures  of  interest.  Minor  corrosion  (<  10X  hole  enlargement)  was  noted  for  BF 

and  H2.  J 

(2)  Under  the  experimental  conditions  discussed,  methane  pyrolyzes 
to  form  large  quantities  of  amorphous  carbon.  The  carbon  particles  were  found 
to  be  spherical  and  ranged  in  size  from  0.1  to  0.05  microns.  No  major  problems 
were  encountered  with  respect  to  clogging  of  0.040  inch  capillaries. 

(3)  Where  the  use  of  segmented  sections  such  as  pyrolytic  graphite 
washers  is  contemplated,  the  alignment  of  the  capillary  holes  is  critical. 


I 
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3.2  COLD  FLOW  TESTS 


a.  Introduction 

There  are  a  number  of  parameters  that  effect  the  efficiency  of  mass 
injection  in  thermal  and  chemical  protection  of  nozzle  walls.  The  interrelation 
between  these  parameters  and  the  complicated  behavior  of  the  fluid  mechanics  of  a 
3-dimensional  wall  jet  make  a  mathematical  treatment  imprac t ica 1 (References  4 
and  37).  In  order  to  investigate  the  dependence  of  important  parameters  on  the 
fluid  dynamic  properties  of  a  wall  jet,  a  cold  flow  experiment  was  designed.  The 
results  obtained  from  these  cold  flow  tests  will  aid  in  the  design  of  an  injection 
hole  geometry  to  obtain  an  efficient  wall  jet. 

From  a  fluid  dynamic  approach,  an  efficient  wall  jet  design  is  one  that 
will  cover  a  large  amount  of  nozzle  wall  area  with  a  protective  film.  With  this 
in  mind,  the  important  parameters  effecting  wall  coverage  of  a  3-dimensional  jet 
(i.e.  jet  spreading  and  downstream  protection)  are:  (1)  hole  spacing,  (2)  hole  size, 
(3)  ,  (axial  location  in  nozzle  of  injection  holes,  (4)  ^U,  >U  ,  and  U  ratios 
of  jet  to  free  stream,  and  (5)  hole  angle.  It  may  be  noted,  that  the  first  four 
of  the  above  parameters  are  all  interdependent  and  difficulty  may  arise  in  their 
separation. 

b.  Apparatus 


The  3-dimensiona 1  conditions  in  the  nozzle  are  simulated  by  the  use  of 
a  1-1/2"  diameter  clear  plastic  tube,  in  which  fully  developed  turbulent  flow  at 
the  point  of  injection  is  insured  by  a  six  foot  entrance  section. 

2 

The  p  U  profiles  are  obtained  with  a  hypodermic  needle  mounted  on  a 
copper  tube.  The  needle  was  ground  to  a  0.011"  I.D.  and  0.014"  O.D.  and  excepted 
only  after  inspection  under  a  microscope.  The  needle  and  copper  tube  (Pitot 
probe)  is  held  in  position  by  two  metal  collars,  which  in  turn  are  supported  by 
four  radial  positioning  studs.  The  positioning  studs  are  then  bolted  to  an  end 
section  for  rigidity.  The  metal  collars  allow  the  Pitot  probe  axial  and  circum¬ 
ferential  degrees  of  freedom.  The  needle's  location  with  respect  to  the  wall  and 
injection  hole  is  determined  from  (1)  radial  pointer  allowing  axial  and  circum¬ 
ferential  measurements  and  (2)  a  dial  indicator  attached  to  one  positioning  stud 
for  radial  distance  from  the  wall. 

The  Pitot  tube^and  a  static  pressure  tap  are  connected  to  a  water  mano¬ 
meter  resulting  in  a  measurement  of  the  flow.  The  water ?manome ter  is  inclined 

17  from  the  horizontal  to  give  an  accuracy  of  +  5  lb/ft-sec  in  p  lr  measurement. 

Injection  gases  with  different  densities  are  used  in  the  experiment 

(1)  study  the  effect  of  density  gradient  on  the  fluid  mechanics  of  a  jet,  and 

(2)  provide  a  method  of  varying  injectant  mass  and  momentum  rates  when  the  velocity 
is  held  constant.  To  provide  a  wide  range  of  injection  densities  the  following 
injection  gases  are  used  in  a  nitrogen  main  stream:  (1)  helium,  (2)  ammonia, 

(3)  nitrogen  and  (4)  freon-22.  It  may  be  noted  here,  that  0 is  chosen  as 
the  fluid  dynamic  parameter  in  comparing  different  experiments.  This  .noice  is 
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based  on  two  reasons:  (1)  the  fact  that  no  concentration  measurements  are  taken 
*n  the  cold  flow  tests,  and  (2)  the  magnification  in  velocity  profile.  That  is, 
fcne  variation  of  the  velocity^ through  the  boundary  layer  is  magnified  due  to  the 
squared  velocity  term  in 

Jet  wall  coverage  can  be  determined  in  one  of  two  ways.  The  first 
method  utilizes  the  (.  u  profiles  near  the  wall;  but  due  to  the  small  magnitudes  in 
and  the  error  in  probe  location,  the  relative  jet  spreading  at  the  wall 
cannot  be  determined  accurately.  A  much  simpler  and  more  accurate  scheme  is  to  line 
the  wall  surface,  downstream  of  injection,  with  Ozalid  paper  and  NH^  as  the  injection 
gas.  This  results  in  a  photographic  method  of  determining  wall  coverage. 

c .  Resu  Its 


In  the  first  quarter  of  the  contract  period,  a  compre!  jnsive  study 
of  a  single  3-dimensional  wall  jet  was  achieved  with  the  cold  flow  apparatus. 

The  study  consisted  of  (1)  the  investigation  of  jet  detachment  (jet  separation 
from  wall),  (2)  development  of  the  jet's  ?  U2  profile  with  axial  and  circumferential 
positions,  and  (3)  dependence  of  hole  geometry  on  wall  coverage.  A  discussion  of 
the  conclusions  obtained  in  the  above  areas  is  presented  below. 


In  the  area  of  jet  detachment,  it  was  required  to  substantiate  the  fact 
that  the  radial  momentum  of  the  wall  jet  in  the  actual  nozzle  did  not  induce  separa¬ 
tion.  This  was  accomplished  by  experimentally  obtaining  maximum  radial  momentum 
rates  that  insured  jet  attachment  for  given  free  stream  velocities.  QFor  a  constant 
main  stream  flow  and  injection  hole  diameter  of  0.0625,f  (inclined  45  from  the  tube 
axis),  the  jet's  momentum  rate  was  varied  and  jet  separation  was  determined  from  a 
p  traverse.  It  may  be  noted  that  the  momentum  rate  of  the  jet  was  varied  in 
one  of  two  ways:  (3)  utilizing  gases  of  different  densities,  and  (2)  changing  the 
injection  mass  flow. 


Typical  0  U  profiles  for  the  attached  and  detached  jet  are  given 
in  Figures  43  and  44  ,  respectively,  with  measurements  taken  at  the  circum¬ 

ferential  center  of  the  jet.  From  the  figure  one  may  conclude  that  near  injection, 
the  p  UZ  profile  for  the  detached  jet,  can  be  divided  into,  essentially,  three 
regions.  The  first  region  is  near  t he  wall  in  which  the  velocity  is  practically 
zero.  The  deveiop.ucrt  of  this  stagnation  region  may  be  compared  to  the  stagnation 
region  at  the  back  side  of  a  cylinder  in  cross  flow;  that  is,  when  detachment  occurs 
the  Jet  (near  the  injection  point)  behaves  as  a  cylinder  in  cross  flow.  The  second 
region  consists  of  the  injected  gas,  which  may  be  called  a  free  jet.  This  free 
jet  region  is  bounded  on  both  sides  by  the  main  stream  gas  and  is  the  "cylinder" 
referenced  previously.  The  third  flow  region  is  that  of  the  main  stream.  Between 
these  three  regions  are  areas  of  large  momentum  transfer  (large  shearing  stress) 
and,  therefore,  the  p  profile  changes  rapidly  in  the  axial  direction.  It  may 
be  of  interest  to  notice  the  redevelopment  of  the  profile  in  the  region  downstream 
from  where  the  stagnation  region  persisted.  This  region  is  strongly  deficient  in 
momentum  and  due  to  the  3-dimentiona 1  shape  of  the  jet,  the  flow  pattern  of  the  main 
stream  around  the  jet  tends  to  develop  a  profile  of  its  own. 


The  attached  wall  jet  profile  is  separated  from  that  of  the  main  stream 
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by  a  free  shear  layer  (Plpjre  43  ).  The  shearing  stress  acting  upon  both  the 

main  stream  and  jet  tends  to  re-establish  the  flow  into  that  of  a  fully  developed 
turbulent  profile.  This  down-stream  re-establishment  region  is  usually  called 
"the  mixing  region"  of  a  wall  Jet. 

In  Table  III  is  presented  the  results  of  jet  detachment,  in  which  the 
radial  momentum  rate  of  the  Jet  is  defined  as  (QU  ),  where  U  is  the  velocity  in 
the  radial  direction  and  Q  is  the  injection  mass  flow  per  unft  area.  Since  the 
velocity  boundary  Ijyer  completely  fills  the  tube  in  developed  pipe  flow,  the 
mementum  rate  (  ?  U  )  of  the  free  stream  was  taken  at  a  height  from  the  wall  equal 
to  the  radius  of  the  injection  hole  (the  average  free  stream  axial  momentum  rate 
acting  upon  the  wall  jet).  From  the  table  a  maximum  radial  momentum  ratio  of  approxi¬ 
mately  0.40  is  obtained  for  jet  attachment.  The  resulting  radial  momentum  rate  is 
applied  to  the  nozzle  only  as  an  order  of  magnitude  in  predicting  jet  separation, 
as  some  questions  may  arise  as  to  the  validity  of  its  direct  application  to  the 
actual  nozzle.  For  example,  from  References  4  and  38,  the  fluid  mechanics  of  a 
wall  jet  is  dependent  on  (1)  temperature  gradient,  (2)  pressure  gradient  (3)  com- 
pressibility  effects,  all  of  which,  were  not  included  in  the  cold  flow  tests. 


TABLE  III 

CORRELATION  OF  JET  DETACHMENT  WITH  MOMENTUM  RATE 


GAS 

Q«r/(  ?u2)^ 

Experimental 

V< 

Position  of  Jet 
with  Respect  to 
wall 

(QUr/pU2)^ 

Theoretical 

N2 

3.60 

2.26 

Detached 

0.60 

He 

3.18 

0.80 

Detached 

0.60 

N2 

3.14 

2.10 

Detached 

0.60 

Freon-22 

1.94 

2.90 

Detached 

0.60 

He 

0.93 

0.44 

Detached 

0.60 

He 

0.38 

0.28 

Attached 

0.60 

Freon-22 

0.3d 

1.28 

Attached 

0.60 

N2 

0.24 

0.53 

Attached 

0.60 
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There  have  been  devised  several  mathematical  solutions  to  the  problem 
of  jet  detachment  in  the  field  of  thrust  vector  control.  An  analysis  which  is 
applicable  to  our  case  is  that  of  Zukowski  and  Spaid,  Reference  39  .  In  their 

solution,  one  assumes  the  shape  of  the  3-dimensional  wall  jet  near  injection  (i.e., 
quarter  sphere  followed  by  an  axisymmetric  half  body)  and  then  calculates  the  pres¬ 
sure  forces  along  the  assumed  body  shape.  Applying  the  conservation  of  momentum 
to  the  radial  direction,  the  radial  pressure  force  equals  the  radial  momentum  rate 
of  the  jet.  With  this  expression,  either  one  of  two  parameters  may  be  found.  The 
first  parameter  is  that  of  jet  detachment.  By  assuming  a  body  shape  with  accom¬ 
panying  dimensions,  jet  detachment  may  be  determined  from  the  fact  that  the  radial 
pressure  force  is  not  great  enough  to  eliminate  the  radial  momentum  rate  of  the  jet. 
The  second  parameter  is  that  of  body  size.  If  the  body  shape  is  assumed,  the  body 
dimensions  may  be  found  by  equating  the  radial  pressure  force  and  momentum  rate. 

This  method  will  be  used  in  determining  the  jet  height  in  the  actual  nozzle. 

To  obtain  a  simple  analysis  of  jet  separation,  similar  to  the  method  of 
Zukoski  and  Spaid,  a  cylinder  in  cross  flow  was  assumed.  Since  the  pressure  dis¬ 
tribution  or  drag  on  a  cylinder  is  very  easily  found  Reference  40  ,  the  radial 
component  of  the  drag  was  equated  to  the  radial  momentum  of  the  injection  gas.  Jet 
detachment  occurs  if  the  radial  momentum  rate  exceeded  the  radial  drag.  This  analy¬ 
sis  is  very  approximate  in  that  cylinder  height  had  to  be  assumed,  which  was  found 
from  experiment,  in  determining  the  drag  force.  The  numerical  results  are 
in  Table  HI  by  the  tabulation  of  the  maxinum  momentum  rate  ratio,  (QU  / 

(see  page >7  for  definition)  that  can  exist  for  jet  attachment.  r 

In  the  study  of  jet  detachment  that  was  carried  out  in  the  first  quarter, 
the  phenomenon  of  downstream  jet  reattachment  was  not  included.  That  is,  due  to 
Jet  expansion  a  detached  jet  may  reattach  to  the  wall  downstream  of  injection. 

This  phenomenon  will  be  studied  in  a  very  simple  manner  in  the  next  quarter. 

2 

The  development  of  the  U  profile  with  both  axial  and  circumferential 
position  was  obtained  in  the  cold  flow  tests  in  order  to  acquire  knowledge  of  the 
complicated  jet  mixing  region  and,  indirectly,  aid  in  a  boundary  layer  analysis.  A 
mathematical  treatment  of  the  mixing  region  results  in  a  set  of  non-linear  differen¬ 
tial  equations.  Reference  37,  with  a  very  costly  and  time  consuming  solution.  There¬ 
fore,  by  performing  a  simple  cold  flow  experiment  the  radial,  axial  and  circumference 
p  distribution  of  a  3-dimension  jet  was  obtained  in  the  mixing  region. 

2 

Radial  and  axial  p  U  profiles  are  presented  in  Figure  45  ,  while 

the  circumferential  traverses  are  jhcwn  in  Figure  46  .  The  injected  gas  was 

Freon-22  with  a  hole  size  of  0.070  and  injected  20°  from  the  tube  axis.  It  may 
be  noted  that  a  similar  experiment  was  performed  using  as  the  injected  gas 

(identical  p>U  ,  hole  size  and  angle)  with  the  f  U2  profiles  of  falling  pre¬ 
cisely  on  those  of  Freon-22.  It  cannot  be  concluded  that  the  ^  U2  profile  is  in¬ 
dependent  of  gas  type  until  further  study  has  been  carried  out.  Although  it  is 
suspected  that,  in  the  mixing  region,  the  flow  field  of  a  wall  jet  is  relatively 
independent  of  gas  type,  since  (1)  diffusion  velocities  are  small,  compared  to 
axjal  velocities,  (2)  difference  in  the  eddy  viscosity  £  *  constant  (dU/dy)  / 

(d  U/dy2)  is  negligible  between  different  gases;  and  (3)  equivalent  turbulent 
mixing  at  injection  occurs  (i.e.  approximately  equal  injection  Reynold's  Number) 


presented 
U  ^theory 
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With  the  effect  of  f>  U  on  velocity  magnification,  discussed  earlier, 
in  mind;  the  mixing  region  may  be  considered  to  terminate  at  an  x/d  between  6 
and  16  (elimination  of  the  inflection  point  in  the  velocity  profiles).  Seban 
and  Korst  Reference  4  obtained,  analytically  and  experimentally,  the  downstream 
location  at  which  the  wall  shear  and  free  jet  boundary  merge;  which  in  essence 
is  only  another  definition  of  the  nixing  region.  When  the  flow  properties  of 
Figure  45  are  applied  to  the  results  of  Seban  and  Korst,  one  obtains  an  x/d 
of  approximately  13  for  the  termination  of  the  mixing  region. 

Probably  the  most  important  investigation  performed  by  the  cold  tests 
is  the  interrelation  of  hole  geometry  on  wall  coverage.  To  obtain  efficient 
chemical  and  thermal  protection  of  a  nozzle  wall,  the  film  must  cover  the  maxi¬ 
mum  amount  of  nozzle  surface.  In  order  to  achieve  this,  the  amount  of  jet  spread¬ 
ing  must  be  established  and  then  the  hole  size,  shape,  spacing  and  axial  lccatinn 
designed  accordingly.  An  analysis  of  jet  spreading  is  no  simple  matter,  as  the 
majority  of  spreading  occurs  in  the  complicated  jet  mixing  region.  Therefore, 
by  lining  the  tube  wall  with  Ozalid  paper  and  using  NH^as  the  injection  gas,  an 

idea  of  wall  coverage  is  obtained  in  isothermal  flow.  This  scheme  is  currently 
in  progress  with  investigation  being  centered  in  areas  of  jet  interaction,  hole 
spacing,  radial  angle  of  injection  (i.e.  angle  producing  a  velocity  in  the  cir¬ 
cumferential  direction  of  a  nozzle)  and  free  stream  boundary  layer  thickness. 

As  noted  before,  direct  application  of  the  cold  flow  results  to  the 
actual  nozzle  is  not  possible.  This  is  due  to  the  fact  that  jet  expansion  is 
very  dependent  on  injection  and  free  stream  £U,  (JU2  and  U  ratios;  Reynold's 
numbers;  density  and  thermal  gradient  between  film  and  free  stream;  momentum 
and  displacement  thickness  of  main  stream  at  injection,  and  free  stream  grad¬ 
ients.  Therefore,  it  is  proposed  that  in  the  next  quarter  of  the  contract  period 
a  study  be  made  to  determine  the  effect  of  the  above  parameters  on  thermal  and 
chemical  protection  of  the  actual  nozzle. 
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SECTION  4 


ROCKET  MOTOR  TESTS 


4.1  CHEMICAL  DIFFUSION  TESTS 


A  technique  has  been  devised  to  study  the  relative  performance  of  in¬ 
jection  fluids.  Basically,  it  consists  of  a  source  of  high-temperature  corrosive 
gas,  a  gas-injection  system,  and  a  cylindrical  test  section.  The  hot-gas  source 
is  a  small  hydrogen-oxygen  gas  generator.  The  gas-injection  system  consists  of  the 
gas  supply,  a  means  of  Internal  distribution,  and  two  injection  slots.  The  in¬ 
jection  fluid  source  will  be  a  high-pressure  gas  bottle  for  fluids  which  are  gases 
at  the  storage  pressures  and  temperatures.  The  cylindrical  test  specimen  is  simply 
a  metallic  tube,  the  inner  surface  of  which  has  been  polished.  The  polished  surface 
will  act  as  a  chemical  reaction  indicator  by  changes  in  its  surface  reflectivity  and 
by  color  changes  caused  by  the  formation  of  the  various  metal  oxides.  During  a  test, 
two  quadrants  of  the  tube  will  be  exposed  directly  to  the  hot  combustion  products 
to  provide  a  check  on  the  injector,  to  observe  the  circumferential  expansion  of  the 
injected  film,  and  to  provide  a  direct  visual  comparison  of  the  effectiveness  of 
the  gas  injections.  The  appearance  of  the  other  two  quadrants  of  the  test  section 
should,  for  stable  film  injection,  vary  from  its  original  appearance  to  something 
close  to  that  of  the  unprotected  sections  at  sufficient  distances  from  the  point 
of  injection.  The  chemical  insulation  distance  can  be  determined  visually  and 
checked  by  microscopic  examination.  The  instrumentation  will  include  measurement 
of  the  generator  chamber  pressure  and  two  sets  of  temperature  measurements  on  the 
outer  surface  of  the  protected  and  unprotected  quadrants  of  the  metal  tube  inserts. 


Most  of  the  tests  will  be  run  under  standard  operating  conditions  which 
are  estimated  to  be  the  following: 


(1) 

Combustion  temperature 

1800°F 

(2) 

Combustion  pressure 

600  psig 

(3) 

Operating  time 

JO -60  sec 

(4) 

Injection  slot  width 

5-10  mils 

(3) 

Tube  length 

6  inches 
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Adjustments  in  any  or  all  of  the  above  may  be  required  to  achieve  reasonable 
chemical  insulation  lengths,  general  economy,  and  good  oxidation  comparisons. 
Deviations  from  the  standard  conditions  wili  be  made  to  evaluate  changes  in 
both  injection  and  combustion  parameters. 

As  mentioned  earlier,  it  is  planned  *  at  the  polished  inner  surface 
of  the  metallic  tube  will  act  as  hemical  reaction  indicator  by  color  changes 
caused  by  the  formation  of  metal  oxide.  Nicke?  >xide,  NiO,  is  black  and  the 
copper  oxides  are  reddish-brown  Cu20,  and  hi  CuO.  In  order  to  st^dy  the 
effectiveness  and  sensitivity  of  nickel  and  '.opper  as  indicators,  several  small 
scale  experiments  were  made  in  which  nickel  and  copper  were  subjected  at  200  psia 
to  small  partial  pressures  of  oxygen.  The  temperature  range  studied  was  600°F 
to  900  F.  The  exposure  times  at  temperature  ranged  from  0.5  min.  to  3  rain,  and 
the  oxygen  concentration  ranged  from  0.5  to  10  volume  pe. cent.  Oxyge.  was  chosen 
rather  than  gaseous  H^O  on  the  basis  of  thermochemical  calculations,  vnich  showed 
H^O  to  be  more  stable^than  either  NiO,  Cu20  or  CuO  at  the  we^eratures  of  interest. 
Tnus,  the  hot  gas  stream  would  have  to  contain  an  excess  of  oxygen  in  order  for 
the  metal  to  act  as  an  indicator. 

Under  the  experimental  conditions  of  temperature,  oxygen  concentration 
and  short  exposure  times,  NiO  did  not  form.  On  the  other  hand,  above  900°F  copper 
oxidized  readily,  even  at  low  oxygen  concentrations,  0.5  volume  ^rcent.  and  short 
exposure  times,  0.5  min.  On  this  basis,  copper  has  been  chosen  a:>  vhe  chemical 
reaction  indicator  for  the  initial  gas-injection  studies.  The  test  nardware  has 
been  completed  and  the  initial  firing  has  been  scheduled  for  the  first  week  of  the 
second  quarter. 


4.2  SIMULATOR  TESTS  (SUB-SCALE) 

a .  Solid  Propellant  Simulator 

The  Simulator  is  a  rocket  motor  which  utilizes  gaseous,  liquid  and  slurried 
propellants  to  simulate  the  exhaust  products  and  temperatures  of  a  given  so? id  pro¬ 
pellant.  In  this  program,  the  requirement  was  to  simulate  a  high  energy  sol  1  con¬ 
taining  16X  aluminum  and  having  a  combustion  temperature  between  6600°F  and  6?'i.)0°F. 

A  chamber  pressure  of  700  psia  was  required. 

0  The  simulation  used  produces  a  theoretical  combustion  temperature  of 
6720  F  with  an  equilibrium  C*  of  5320  fps.  Assuming  97X  combustion  efficiency  a 
total  flow  rate  of  5.35  lb/sec  will  produce  700  psia  chamber  pressure  with  a  1.25  in. 
diameter  throat. 

The  propellant  consists  of  gaseous  oxygen  and  nitrogen  and  a  slurried  fuel 
containing  powdered  aluminum.  Exhaust  products  are  given  in  Table  I  ,  Section  2.1. 

b.  Description  of  First  Sub-scale  Nozzle 

The  first  sub-scale  nozzle  (See  Figure  47  )  will  be  a  pyrolytic  graph¬ 

ite  heat  sink  convectively  cooled  with  methane.  There  will  be  no  film  injection 
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FIGURE  47.  NOZZLE  FOR  FIRS7  SUB-SCALE  TEST 


in  this  test.  The  cooling  passages  will  be  nominally  G.25  inches  from  tn<?  sur¬ 
face  and  follow  the  contour  as  nearly  as  possible.  The  threat  washer  contains 
28-0.030  in.  diameter  holes.  The  upstream  and  downstream  washers  contain  23- 
0.040  in.  diameter  holes.  The  coolant  will  be  introduced  downstream  and  flow  to¬ 
ward  the  nozzle  inlet. 

The  nozzle  throat  diameter  is  1,25  inches  and  the  heat  sink  outside 
diameter  is  5.00  inches. 

The  methane  coolant  flow  rate  of  0.001  Ib/sec/hole  giving  a  total  of 
0.028  lb/sec.  According  to  the  computer  results  the  surface  temperature  at  the 
throat  should  not  exceed  5020  F  in  50  seconds  (the  planned  duration). 

Instrumentation  will  include  chamber  pressure,  propellant  flow  rates, 
thrust,  methane  pressures  to  determine  flow  rate  and  pressure  drops,  methane 
ou^et  temperature,  and  heat  sink  temperatures  at  five  selected  locations. 
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